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Abstract. 

A review of the experimental and theoretical determinations of the anomalous 
magnetic moment of the muon is given. The anomaly is defined by a = (g — 2)/2, where 
the Lande g-factor is the proportionality constant that relates the spin to the magnetic 
moment. For the muon, as well as for the electron and tauon, the anomaly a differs 
slightly from zero (of order 10~^) because of radiative corrections. In the Standard 
Model; contributions to the anomaly come from virtual 'loops' containing photons and 
the known massive particles. The relative contribution from heavy particles scales as 
the square of the lepton mass over the heavy mass, leading to small differences in the 
anomaly for e, and r. If there are heavy new particles outside the Standard Model 
which couple to photons and/or leptons, the relative effect on the muon anomaly will 
be ^ {rai^JmeY' w 43 x 10'^ larger compared with the electron anomaly. Because 
both the theoretical and experimental values of the muon anomaly are determined 
to high precision, it is an excellent place to search for the effects of new physics, or 
to constrain speculative extensions to the Standard Model. Details of the current 
theoretical evaluation, and of the series of experiments that culminates with E821 at 
the Brookhaven National Laboratory are given. At present the theoretical and the 
experimental values are known with a similar relative precision of 0.5 ppm. There is, 
however, a 3.4 standard deviation difference between the two, strongly suggesting the 
need for continued experimental and theoretical study. 
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1. Introduction and History of (7-Factors 

A charged elementary particle with half-integer intrinsic spin has a magnetic dipole 
moment /I aligned with its spin s: 



where q = ±e is the charge of the particle in terms of the magnitude of the electron 
charge e, and the proportionality constant Qs is the Lande (yf-factor. For the charged 
leptons, e, /i, or r, Qg is slightly greater than 2. From a quantum-mechanical view, the 
magnetic moment must be directed along the spin, since in the absence of an external 
electric or magnetic field, the spin provides the only preferred direction in space. 

While the behavior of the spin of an individual particle must be treated quantum 
mechanically, the polarization, the average behavior of the spins of a large ensemble 
of particles, can be treated to a large extent as a classical collection of spinning bar 
magnets. II An externally applied magnetic field, B, will exert a torque on the magnetic 
moment which tends to align the polarization with the direction of the field. However, 
because the magnetic moments are associated with the intrinsic angular momentum, the 
polarization precesses in this case like a classical collection of gyroscopes in the plane 
perpendicular to the field. In both the classical and quantum-mechanical cases, the 
torque on the particle is given hj N = fl x B, the energy depends on the orientation of 
the magnetic dipole, Ti. = —fl ■ B, and there is a net force on the particle if the field is 
non-uniform. 

The study of atomic and sub-atomic magnetic moments, which continues to this 
day, began in 1921 with the famous " Stern- Gerlach" experiments [1, 2]. A beam of 
silver atoms was passed through a gradient magnetic field to separate the individual 
quantum states. Two spatially separated bands of atoms were observed, signifying two 
quantum states. From this separation, the magnetic moment of the silver atom was 
determined to be one Bohr magneton[3], eh/2me, to within 10%. Phipps and Taylor[4] 
repeated the experiment with a hydrogen beam in 1927, and they also observed two 
bands. From the splitting of the bands they concluded that, like silver, the magnetic 
moment of the hydrogen atom was one Bohr magneton. In terms of spin, which was 
proposed independently by Compton[5], and by Uhlenbeck and Goudsmit[6], the two- 
band structure indicated that the ^-component of the angular momentum had two 
values, Sz = ±h/2. Later, the magnetic moments of both atoms would be traced to the 
un-paired spin of an atomic electron, implying (from Equation 1) that the (7-factor of 
the electron is 2. 

In 1928 Dirac published his relativistic wave equation for the electron. It employed 
four-vectors and the famous 4x4 matrix formulation, with the spin degree of freedom 
emerging in a natural way. Dirac pointed out that his wave equation for an electron in 

11 Following the custom in the literature, in this article we will often refer to the motion of the 'spin' 
when more correctly we are speaking of the motion of the polarization of an ensemble of a large number 
of particles which have spin. 




(1) 
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external electric and magnetic fields has "the two extra terms^ 

-(cr,H)+z-pi(c^,E), (2) 

c c 

. . . when divided by the factor 2m can be regarded as the additional potential energy 
of the electron due to its new degree of freedom[7]." These terms represent the 
magnetic dipole (Dirac) moment and electric dipole moment interactions with the 
external magnetic and electric fields. Dirac theory predicts that the electron magnetic 
moment is one Bohr-magneton [qs = 2), consistent with the value measured by the 
Stern-Gerlach and Phipps- Taylor experiments. Dirac later commented: "It gave just 
the properties that one needed for an electron. That was an unexpected bonus for me, 
completely unexpected[8]." The non-relativistic reduction of the Dirac equation for the 
electron in a weak magnetic field [9], 

d(f) 



ih-^ 

dt 



'L + 2S]-B 



(3) 



2m 2m 

shows clearly that Qs = 2, and the gf-factor for orbital angular momentum, ge = 1. 

With the success of the Dirac equation, it was believed that the proton should also 
have a (^-factor of 2. However, in 1933 Stern and his collaborators [10] showed that the 
(yf-factor of the proton was ~ 5.5. In 1940 Alvarez and Bloch[ll] found that the neutron 
likewise had a large magnetic moment, which was a complete surprise, since g = for 
the neutron. These two results remained quite mysterious for many years, but we now 
understand that the baryon magnetic moments are related to their internal quark-gluon 
structure. 

Theoretically it is useful to break the magnetic moment into two pieces: 

u = (1 + a) — — where a = . (4) 

^ ^ ^ 2m 2 ^ ^ 

The first piece, predicted by the Dirac equation and called the Dirac moment, is 1 in 

units of the appropriate magneton, qh/2m. The second piece is the anomalous (Pauli) 

moment [12], where the dimensionless quantity a is referred to as the anomaly. 

In 1947, motivated by measurements of the hyperfine structure in hydrogen that 

obtained splittings larger than expected from the Dirac theory [13, 14, 15], Schwinger[16] 

showed that from a theoretical viewpoint these "discrepancies can be accounted for by 

a small additional electron spin magnetic moment" that arises from the lowest-order 

radiative correction to the Dirac moment.^ This radiative correction, which we now call 

the one-loop correction to = 2, is shown diagrammatically in Figure 1(b). Schwinger 

obtained the value Oe = a/(2-7r) ^ 0.00116- ■■ (which is also true for and Qt)- In 

the same year, Kusch and Foley [18] measured with 4% precision, and found that the 

measured electron anomaly agreed well with Schwinger's prediction. In the intervening 

^ This expression uses Dirac's original notation. 

+ In response to Nafe, et al.[13], Breit[17] conjectured that this discrepancy could be explained by the 
presence of a small Pauli moment. It's not clear whether this paper influenced Schwinger's work, but in 
a footnote Schwinger states: "However, Breit has not correctly drawn the consequences of his empirical 
hypothesis." 
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time since the Kusch and Foley paper, many improvements have been made in the 
precision of the electron anomaly [19]. Most recently has been measured to a relative 
precision of 0.65 ppb (parts per billion) [20], a factor of 6 improvement over the famous 
experiments at the University of Washington [21]. 




(a) (b) (c) 

Figure 1. The Feynman graphs for: (a) g = 2; (b) the lowest-order radiative correction 
first calculated by Schwinger; and (c) the vacuum polarization contribution, which is 
an example of the next-order term. The * emphasizes that in the loop the muon is 
off-shell. 

The ability to calculate loop diagrams such as those shown in Figure 1 (b) and (c) 
is intimately tied to the renormalizability of the theory, which provides a prescription 
to deal with the infinities encountered in calculating radiative corrections, and was 
important to the development of quantum electrodynamics. In his original paper 
Schwinger [16] described a new procedure that transformed the Dirac Hamiltonian to 
include the electron self-energy which arises from the emission and absorption of virtual 
photons. By doing so, he eliminated the divergences encountered in calculating the 
lowest-order radiative correction. He pointed to three important features of his new 
Hamiltonian: "it involves the experimental electron mass" (known today as the 'dressed 
or physical mass'); "an electron now interacts with the radiation field only in the presence 
of an external field;" (i.e. the virtual photons from the self-energy are absent) and "the 
interaction energy of an electron with an external field is now subject to a finite radiative 
correction" [16]. This concept of renormalization also played an important role in the 
development of the Standard Model, and the lowest-order contribution from virtual W 
and Z gauge bosons to was calculated very soon after the electroweak theory was 
shown to be renormalizable[22]. 

The diagram in Figure 1(a) corresponds to (7 = 2, and the first-order (Schwinger) 
correction which dominates the anomaly is given in diagram 1(b). More generally, 
the Standard-Model value of the electron, muon or tauon anomaly, a(SM), arises from 
loops (radiative corrections) containing virtual leptons, hadrons and gauge bosons. By 
convention, these contributions are divided into three classes: the dominant QED terms, 
like Schwinger's correction, which contain only leptons and photons; terms which involve 
hadrons, particularly the hadronic vacuum polarization correction to the Schwinger 
term; and electroweak terms, which contain the Higgs, W and Z. Some of the terms 
are identical for all three leptons, but as noted below, there are mass-dependent terms 
which are significant for the muon and tauon but not for the electron. As a result, the 
muon anomaly is slightly larger than that of the electron. 
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Both the QED and electroweak contributions can be calculated to high precision. 
For the muon, the former has been calculated to 8th order (4th order in a/vr) and the 
10th order terms are currently being evaluated (see §7.1). The electroweak contributions 
have been calculated to two loops, as discussed in §7.4. 

In contrast, the hadronic contribution to cannot be accurately calculated from 
low-energy quantum chromodynamics (QCD), and contributes the dominant theoretical 
uncertainty on the Standard-Model prediction. The lowest-order hadronic contribution, 
determined using experimental data via a dispersion relation as discussed in §7.2, 
accounts for over half of this uncertainty. The Hadronic light by light scattering 
contribution (see §7.3) contributes the rest. This latter contribution cannot be related to 
existing data, but rather must be calculated using models that incorporate the features 
of QCD. 

Since the value of arises from all particles that couple, directly or indirectly, to 
the muon, a precision measurement of the muon anomaly serves as an excellent probe for 
new physics. Depending on their mass and coupling strength, as yet unknown particles 
could make a significant contribution to a^. Lepton or IV-boson substructure might 
also have a measurable effect. Conversely, the comparison between the experimental 
and Standard-Model values of the anomaly can be used to constrain the parameters of 
speculative theories [23, 24, 25]. 

While the current experimental uncertainty of ±0.5 ppm (parts per million) on the 
muon anomaly is 770 times larger than that on the electron anomaly, the former is 
far more sensitive to the effects of high mass scales. In the lowest-order diagram where 
mass effects appear, the contribution of heavy virtual particles scales as (miepton/''^Hv)^, 
giving the muon a factor of (m^/me)^ ~ 43000 increase in sensitivity over the electron. 
Thus, at a precision of 0.5 ppm, the muon anomaly is sensitive to physics at a few 
hundred GeV scale, while the reach of the electron anomaly at the current experimental 
uncertainty of 0.65 ppb is limited to the few hundred MeV range. 

To observe effects of new physics in Oe, it is necessary to have a sufficiently 
precise standard-model value to compare with the experimental value. Since the QED 
contribution depends directly on a power series expansion in the fine-structure constant 
a, a meaningful comparison requires that a be known from an independent experiment 
to the same relative precision as ag. At present, independent measurements of a have 
a precision of ~ 6.7 ppb [27], about ten times larger than the error on Og. On the other 
hand, if one assumes that there are no new physics contributions to Cg, i.e. the Standard- 
Model radiative corrections are the only source of Cg, one can use the measured value 
of Og and the standard-model calculation to obtain the most precise determination of 
a [28]. 

1.1. Muon Decay 

Since the kinematics of muon decay are central to the measurements of a^, we discuss 
the general features in this section. Specific issues that relate to the design of E821 at 
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Brookhaven are discussed in the detector section, § 2.7.1. 

After the discovery of parity violation in /3-decay[29] and in muon decay[30, 31], it 
was reahzed that one could make beams of polarized muons in the pion decay reactions 

TT" ^ /i" + or 7r+ ^ /i+ + v^. 

The pion has spin zero, the neutrino (antineutrino) has helicity of -1 (+1), and the 
forces in the decay process are very short-range, so the orbital angular momentum in 
the final state is zero. Thus conservation of angular momentum requires that the /i~ 
(/X"'") helicity be +1 (-1) in the pion rest frame. The muons from pion decay at rest are 
always polarized. 

From a beam of pions traversing a straight beam-channel consisting of focusing 
and defocusing elements (FODO), a beam of polarized, high energy muons can be 
produced by selecting the "forward" or "backward" decays. The forward muons are 
those produced, in the pion rest frame, nearly parallel to the pion laboratory momentum 
and are the decay muons with the highest laboratory momenta. The backward 
muons are those produced nearly anti-parallel to the pion momentum and have the 
lowest laboratory momenta. The forward /i~ (/i^) are polarized along (opposite) their 
lab momenta respectively; the polarization reverses for backward muons. The E821 
experiment uses forward muons, produced by a pion beam with an average momentum 
of ~ 3.15 GeV/c. Under a Lorentz transformation from the pion rest frame to the 
laboratory frame, the decay muons have momenta in the range < < 3.15 GeV/c. 
After momentum selection, the average momentum of muons stored in the ring is 3.094 
GeV/c, and the average polarization is in excess of 95%. 

The pure {V — A) three-body weak decay of the muon, yU~ ^ + z/^ + or 
/i e+ -|- z/^ -|- z/g, is "self-analyzing", that is, the parity-violating correlation between 
the directions in the muon rest frame (MRF) of the decay electron and the muon spin 
can provide information on the muon spin orientation at the time of the decay. (In the 
following text, we use 'electron' generically for either e~ and e+ from the decay of the 
.) Consider the case when the decay electron has the maximum allowed energy in 
the MRF, -E^ax ~ ('^AtC^)/2 = 53 MeV. The neutrino and anti-neutrino are directed 
parallel to each other and at 180° relative to the electron direction. The vv pair carry 
zero total angular momentum, since the neutrino is left-handed and the anti-neutrino is 
right-handed; the electron carries the muon's angular momentum of 1/2. The electron, 
being a lepton, is preferentially emitted left-handed in a weak decay, and thus has a 
larger probability to be emitted with its momentum anti-parallel rather than parallel to 
the yU.~ spin. By the same line of reasoning, in decay, the highest-energy positrons 
are emitted parallel to the muon spin in the MRF. 

In the other extreme, when the electron kinetic energy is zero in the MRF, 
the neutrino and anti-neutrino are emitted back-to-back and carry a total angular 
momentum of one. In this case, the electron spin is directed opposite to the muon spin in 
order to conserve angular momentum. Again, the electron is preferentially emitted with 
helicity -1, however in this case its momentum will be preferentially directed parallel to 
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the yU~ spin. The positron, in /i"*" decay, is preferentially emitted with helicity +1, and 
therefore its momentum will be preferentially directed anti-parallel to the spin. 

With the approximation that the energy of the decay electron E' » rrieC^, the 
differential decay distribution in the muon rest frame is given by[32], 

dP{y', e') oc n\y') [1 ± A{y') cos 6'] dy'dO! (5) 

where y' is the momentum fraction of the electron, y' = p'Jpl^ ^nax' is the solid angle, 
B' = cos~^ (pg ■ s) is the angle between the muon spin and p g, p'^ ^^^c ^ E'^^^, and 
the (— ) sign is for negative muon decay. The number distribution n{y') and the decay 
asymmetry A{y') are given by 

niy') = 2y"i3-2y') and Aiy') = f^. (6) 

Note that both the number and asymmetry reach their maxima at y' = 1, and the 
asymmetry changes sign at = i, as shown in Figure 2(a). 




Energy, MeV Energy, GeV 

(a) Muon Rest Frame (b) Laboratory Frame 



Figure 2. Number of decay electrons per unit energy, N (arbitrary units), value of 
the asymmetry A, and relative figure of merit NA^ (arbitrary units) as a function of 
electron energy. Detector acceptance has not been incorporated, and the polarization is 
unity. For the third CERN experiment and E821, Emax « 3.1 GeV (p^ = 3.094 GeV/c) 
in the laboratory frame. 

The CERN and Brookhaven based muon [g — 2) experiments stored relativistic 
muons in a uniform magnetic field, which resulted in the muon spin precessing with 
constant frequency Ua, while the muons traveled in circular orbits. If all decay electrons 
were counted, the number detected as a function of time would be a pure exponential; 
therefore we seek cuts on the laboratory observables to select subsets of decay electrons 
whose numbers oscillate at the precession frequency. Recalling that the number of decay 
electrons in the MRF varies with the angle between the electron and spin directions, 
the electrons in the subset should have a preferred direction in the MRF when weighted 
according to their asymmetry as given in Equation 5. At ~ 3.094 GeV/c the 
directions of the electrons resulting from muon decay in the laboratory frame are very 
nearly parallel to the muon momentum regardless of their energy or direction in the 
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MRF. Therefore the only practical remaining cut is on the electron's laboratory energy. 
Typically, selecting an energy subset will have the desired effect: there will be a net 
component of electron MRF momentum either parallel or antiparallel to the laboratory 
muon direction. For example, suppose that we only count electrons with the highest 
laboratory energy, around 3.1 GeV. Let z indicate the direction of the muon laboratory 
momentum. The highest-energy electrons in the laboratory are those near the maximum 
MRF energy of 53 MeV, and with MRF directions nearly parallel to z. There are more 
of these high-energy electrons when the yU~ spins are in the direction opposite to z 
than when the spins are parallel to z. Thus the number of decay electrons reaches a 
maximum when the muon spin direction is opposite to f , and a minimum when they are 
parallel. As the spin precesses the number of high-energy electrons will oscillate with 
frequency uja- More generally, at laboratory energies above ~ 1.2 GeV, the electrons 
have a preferred average MRF direction parallel to z (see Figure 2). In this discussion, 
it is assumed that the spin precession vector, cJ^, is independent of time, and therefore 
the angle between the spin component in the orbit plane and the muon momentum 
direction is given by ujat + 0, where is a constant. 

Equations 5 and 6 can be transformed to the laboratory frame to give the electron 
number oscillation with time as a function of electron energy, 

N,{t, E) = N,o{E)e-"^^[l + cos(a;,t + (7) 

or, taking all electrons above threshold energy Eth, 

N{t, Eth) = N,{Eth)e-"^^[l + A{Eth) cos(a;,t + (8) 

In Equation 7 the differential quantities are, 

ME)=V ~^jf^y^^l , iV,o(i?)oc(2/-l)(V-52/-5), (9) 

and in Equation 8, 

yth{'2yth + 1) 
-Vth + Vth 

In the above equations, y = E/Emax, Vth = Eth/Emax, V is the polarization of 
the muon beam, and E, Eth-, and -Emax = 3.1 GeV are the electron laboratory energy, 
threshold energy, and maximum energy, respectively. 

The fractional statistical error on the precession frequency, when fitting data 
collected over many muon lifetimes to the five-parameter function (Equation 8), is given 
by 

5e=^-^= ^ . . (11) 

where N is the total number of electrons, and A is the asymmetry, in the given data 
sample. For a fixed magnetic field and muon momentum, the statistical figure of merit 
is NA"^, the quantity to be maximized in order to minimize the statistical uncertainty. 



N{E,,) oc {y,, - lf{-yl + y^n + 3), A{E,n) = V ''^ (10) 
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(a) No detector acceptance or energy resolu- (b) Detector acceptance and energy resolution 
tion included included 



Figure 3. The integral A'', A, and N (arbitrary units) for a single energy- 
threshold as a function of the threshold energy; (a) in the laboratory frame, not 
including and (b) including the effects of detector acceptance and energy resolution 
for the E821 calorimeters discussed below. For the third CERN experiment and E821, 



3.1 GcV (p^j = 3.094 GcV/c) in the laboratory frame. 



The energy dependences of the numbers and asymmetries used in Equations 7 and 
8, along with the figures of merit N ^ are plotted in Figures 2 and 3 for the case of 
E821. The statistical power is greatest for electrons at 2.6 GeV (Figure 2). When a fit 
is made to all electrons above some energy threshold, the optimal threshold energy is 
about 1.7-1.8 GeV (Figure 3). 

1.2. History of the Muon {g — 2) Experiments 

In 1957, at the Columbia-Nevis cyclotron, the spin rotation of a muon in a magnetic 
field was observed for the first time. The torque exerted by the magnetic field on the 
muon's magnetic moment produces a spin precession frequency 

qgB qB 

us = -- 1-7, 12 

2m 7m 

where 7 = (1 — with j3 = v/c. Garwin et al.[30] found that the observed rate of 

spin rotation was consistent with g = 2. 

In a subsequent paper, Garwin et al. [33] reported the results of a second experiment, 
a measurement of the muon anomaly to a relative precision of 6.6%, g^ > 2(1.00122 ± 
0.00008), with the inequality coming from the poor knowledge of the muon mass. In a 
note added in proof, the authors reported that a new measurement of the muon mass 
permitted them to conclude that g^ = 2(1.00113to;oooi2)- Within the experimental 
uncertainty, the muon's anomaly was equal to that of the electron. More generally, the 
muon was shown to behave like a heavy electron, a spin 1/2 fermion obeying QED. 

In 1961 the first of three experiments to be carried out at CERN reported a more 
precise result obtained at the CERN synchrocyclotron [34]. In this experiment, highly 
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polarized muons of momentum 90 MeV/c were injected into a 6-meter long magnet 
with a graded magnetic field. As the muons moved in almost circular orbits which 
drifted transverse to the gradient, their spin vectors precessed with respect to their 
momenta. The rate of spin precession is readily calculated. Assuming that f3 ■ B = 0, 
the momentum vector of a muon undergoing cyclotron motion rotates with frequency 

= (13) 

The spin precession relative to the momentum occurs at the difference frequency, Ua, 
between the spin frequency in Equation 12 and the cyclotron frequency, 

- - - f9-2\qB qB 

\ 2 / m m 

The precession frequency Ua has the important property that it is independent of the 
muon momentum. When the muons reached the end of the magnet, they were extracted 
and their polarizations measured. The polarization measurement exploited the self- 
analyzing property of the muon: more electrons are emitted opposite than along the 
muon spin. For an ensemble of muons, Ua is the average observed frequency, and B is 
the average magnetic field obtained by folding the muon distribution with the magnetic 
field map. 

The result from the first CERN experiment was[34] = 0.001 145(22) (1.9%), 
which can be compared with a/2'K = 0.001 161 410 ■■ ■. With additional data this 
technique resulted in the first observation of the effects of the (a/vr)^ term in the QED 
expansion [35]. 

The second CERN experiment used a muon storage ring operating at 1.28 GeV/c. 
Vertical focusing was achieved with magnetic gradients in the storage-ring field. While 
the use of magnetic gradients to focus a charged particle beam is quite common, it makes 
a precision determination of the (average) magnetic field which enters into Equation 14 
rather difficult for two reasons. Since the field is not uniform, information on where the 
muons are in the storage ring is needed to correct the average field for the gradients 
encountered. Also, the presence of gradient magnetic fields broadens the NMR line- 
shape, which reduces the precision on the NMR measurement of the magnetic field. 

A temporally narrow bunch of 10^^ protons at 10.5 GeV/c from the CERN proton 
synchrotron (PS) struck a target inside the storage ring, producing pions, a few of which 
decay in such a way that their daughter muons are stored in the ring. A huge flux of 
other hadrons was also produced, which presented a challenge to the decay electron 
detection system. The electron detectors could only be placed in positions around the 
ring well-removed from the production target, which limited their geometric coverage. 
Of the pions which circulated in the ring for several turns and then decayed, only one in 
a thousand produced a stored muon, resulting in about 100 stored muons per injected 
proton bunch. The polarization of the stored muons was 26% [36]. 

In all of the experiments discussed in this review, the magnetic field was measured 
by observing the Larmor frequency of stationary protons, Up, in nuclear magnetic 
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resonance (NMR) probes. Ignoring the signs of the muon and proton charges, we have 
the two equations: 

ooa = —a^B and ujp = gp \ - — , (15) 
where lo^ is the Larmor frequency for a jree proton. Dividing and solving for we find 

O/. = T = T ^, 16 

A — UJa/UJp A — l-i 

where the fundamental constant A = fifj_/fip, and TZ = Ua/^p- The tilde on uja indicates 
that the measured frequency has been adjusted for any necessary (small) corrections, 
such as the pitch and radial electric field corrections discussed in §3. Equation 16 
was used in the second CERN experiment, and in all of the subsequent muon {g — 2) 
experiments, to determine the value of the anomaly. The most accurate value of the 
ratio A is derived from measurements of the muonium (the fi~^e~ atom) hyperfine 
structure[39, 40], i.e., from measurements made with the //"'", so that it is more properly 
written A+ = /i^j+//ip. Application of Equation 16 to determine a^- requires the 
assumption of CPT invariance. 

The second CERN experiment obtained a^, = (11661.6 ± 3.1) x 10"^ (0.27%) [36], 
testing quantum electrodynamics for the muon to the three-loop level. Initially this 
result was 1.8 standard deviations above the QED value, which stimulated a new 
calculation of the QED light-by-light scattering contribution[lll] that brought theory 
and experiment into agreement. 

In the third experiment at CERN, several of the undesirable features of the second 
experiment were eliminated. A secondary pion beam produced by the PS was brought 
into the storage ring, which was placed a suitable distance away from the production 
target. The vr — * /i decay was used to kick muons onto a stable orbit, producing 
a stored muon beam with polarization of 95%. The "decay kick" had an injection 
efficiency of 125 ppm, with many of the remaining pions striking objects in the storage 
ring and causing background in the detectors after injection. While still significant, 
the background in the electron calorimeters associated with beam injection was much 
reduced compared to the second CERN experiment. 

Another very important improvement in the third CERN experiment, adopted by 
the Brookhaven based experiment as well, was the use of electrostatic focusing[37], which 
represented a major conceptual breakthrough. The magnetic gradients used by the 
second CERN experiment made the precise determination of the average magnetic field 
difficult, as discussed above. The more uniform field greatly improves the performance 
of NMR techniques, to provide a measure of the field at the 10"^ to 10"^ level. With 
electrostatic focusing Equations 12 and 13 must be modified. The cyclotron frequency 
becomes 

q \b 7 (PxE]] 

uJc = 2 T ' ^ 

m7 7"^ — 1\ c / 
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and the spin precession frequency becomes [42] 

r 



m 



1 



B 



1 



7 



/3xE 



7/ /7 + i' ■ yz 7 + ly\ c 

Substituting for = {g^ — 2)/2, we find that the spin difference frequency is 



^a = 

m 



7 



.7 + 1 



If /? ■ -B = 0, this reduces to 



m 



a^B 



iP ■ B)(3 



px E 



px E 



y - 1^ 



1 



(19) 



(20) 



For 7magic = 29.3 (p^ = 3.09 GeV/c), the second term vanishes; one is left with the 
simpler result of Equation 14, and the electric field does not contribute to the spin 
precession relative to the momentum. The great experimental advantage of using the 
magic 7 is clear. By using electrostatic focusing, a gradient B field is not needed for 
focusing and B can be made as uniform as possible. The spin precession is determined 
almost completely by Equation 14, which is independent of muon momentum; all muons 
precess at the same rate. Because of the high uniformity of the B-field, a precision 
knowledge of the stored beam trajectories in the storage region is not required. 

Since the spin precession period of 4.4 fis is much longer than the cyclotron period 
of 149 ns, during a single precession period a muon samples the magnetic field over 
the entire azimuth 29 times. Thus the important quantity to be made uniform is the 
magnetic field averaged over azimuth. The CERN magnet was shimmed to an average 
azimuthal uniformity of ±10 ppm, with the total systematic error from all issues related 
to the magnetic field of ±1.5 ppm[38]. 

Two small corrections to Ua are required to form cDq, which is used in turn to 
determine in Equation 16. 1) The vertical pitching motion (betatron oscillations) 
about the mid-plane of the storage region means that P ■ B differs slightly from zero (see 
Equation refeq:Ediffreq). 2) Only muons with the central radius of 7112 mm are at the 
magic value of 7, so that a radial electric field will effect the spin precession of muons 
with 7 7^ 'jmagic- The Small corrections are described in §3.2. 

The design and shimming techniques of the CERN III storage ring were chosen to 
produce a very stable magnet mechanically. After cycling the magnet power, the field 
was reproducible to a few ppm[56]. The 14 m-diameter storage ring ( Bq = 1.5 T) 
was constructed using 40 identical C-magnets, bolted together in a regular polygon to 
make a ring, and the field was excited by four concentric coils connected in series. The 
coils that excited the field in the CERN storage ring were conventional warm coils, and 
it took about four days after powering for thermal and mechanical equilibrium to be 
reached[56, 57]. The shimming procedure involved grinding steel from the pole pieces, 
which introduced a periodic shape to the contour map averaged over azimuth (see Figure 
5 of Reference [38]). The magnetic field was mapped once or twice per running period, 
by removing the vacuum chamber and stepping NMR probes through the storage region. 
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The electrostatic quadrupoles had a two-fold symmetry, covered approximately 80% 
of the azimuth, and defined a rectangular storage region 120 mm (radial) by 80 mm 
(vertical). The ring behaved as a weak focusing storage ring (betatron) [43, 44, 45], with 
a field index of ~ 0.135 (see §3). 

The pion beam was brought into the storage ring through a pulsed co-axial 
"infiector", which produced a magnetic field that canceled the 1.5 T storage- ring field 
and permitted the beam to arrive undefiected, tangent to the storage circle but displaced 
76 mm radially outward from the central orbit. The infiector current pulse rose to a peak 
current of 300 kA in 12 /is. The transient effects from this pulsed device on the magnetic 
field seen by the stored muons during the precession measurement were estimated to be 
small [38]. 

The decay electrons were detected by 24 lead-scintillator sandwiches each viewed 
through an air light-guide by a single 5-inch photomultiplier tube. The photomultiplier 
signal from each detector was discriminated at several analog thresholds, thus providing 
a coarse pulse height measurement as well as the time of the pulse. Each of the 
discriminated signals was assigned to a time bin, in a manner independent of both the 
ambient rate and the electron arrival time relative to the clock time-boundary. The final 
precision of 7.3 ppm (statistics dominated) convincingly confirmed the effect of hadronic 
vacuum polarization, which was predicted to contribute to at the level of 60 ppm. 
The factor of 35 increase in sensitivity between the second and third experiments can 
be traced to the improved statistical power provided by the injected pion beam, both 
in the number and asymmetry of the muon sample, and in reduced systematic errors 
resulting from the use of electrostatic focusing with the central momentum set to ■jmagic- 

Around 1984 a new collaboration (E821) was formed with the aim of improving 
on the CERN result to a relative precision of 0.35 ppm, a goal chosen because it was 
20% of the predicted first-order electroweak contribution to a^. During the subsequent 
twenty-year period, over 100 scientists and engineers contributed to E821, which 
was built and performed at the Brookhaven National Laboratory (BNL) Alternating 
Gradient Synchrotron (AGS). A measurement at this level of precision would test the 
electroweak renormalization procedure for the Standard Model and serve as a very 
sensitive constraint on new physics[23, 24, 25]. 

The E821 experiment at Brookhaven [47, 48, 49, 50, 51, 26] has advanced the relative 
experimental precision of to 0.54 ppm, with a several standard-deviation difference 
from the prediction of the Standard Model. Assuming CPT symmetry, viz. = a^-, 
the new world average is 

a^(Expt) = 11659 208.0(6.3) X 10"^° (0.54ppm). (21) 

The total uncertainty includes a 0.46 ppm statistical uncertainty and a 0.28 ppm 
systematic uncertainty, combined in quadrature. 

A summary showing the physics reach of each of the successive muon (g-2) 
experiments is given in Table 1. In the subsequent sections we review the most recent 
experiment, E821 at Brookhaven, and the Standard-Model theory with which it is 
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Table 1. Measurements of the muon anomalous magnetic moment. When the 
uncertainty on the measurement is the size of the next term in the QED expansion, 
or the hadronic or weak contributions, the term is hsted under "sensitivity" . The 
"?" indicates a result that differs by greater than two standard deviations with the 
Standard Model. For completeness, we include the experiment of Henry, et al.,[46], 
which is not discussed in the text. 



± 


Measurement 




Sensitivity 


Reference 




g = 2.00 ±0.10 




g = 2 


Garwin et al[30], Nevis (1957) 




001 13+0 00016 

U.UUi iO_Q QQQ]^2 


12.4% 


a 

TV 


Garwin et al[33], Nevis (1959) 




0.001 145(22) 


1.9% 


a 


Charpak et 0^34] CERN 1 (SC) (1961) 




0.001162(5) 


0.43% 




Charpak et al[35] CERN 1 (SC) (1962) 




0.00116616(31) 


265 ppm 




Bailey et al[36] CERN 2 (PS) (1968) 




0.001060(67) 


5.8% 


a 

TT 


Henryei aZ[46] solenoid (1969) 


/^^ 


0.001165 895(27) 


23 ppm 


(^) + Hadronic 


Bailey et al[37] CERN 3 (PS) (1975) 


/^^ 


0.001165911(11) 


7.3 ppm 


(f f + Hadronic 


Bailey et a/[38] CERN 3 (PS) (1979) 




0.001165 9191(59) 


5 ppm 


[^y^ + Hadronic 


Brown et al[48] BNL (2000) 




0.001165 920 2(16) 


1.3 ppm 


(°)4 + Weak 
(f )^ + Weak + ? 


Brown et al[49] BNL (2001) 




0.001165 920 3(8) 


0.7 ppm 


Bennett et al[50] BNL (2002) 




0.001 165 921 4(8)(3) 


0.7 ppm 


(1)4 + Weak + ? 


Bennett et al[5l] BNL (2004) 




0.001165 920 80(63) 


0.54 ppm 


(1)4 + Weak + ? 


Bennett et a/[51, 26] BNL WA (2004) 



compared. 

2. The Brookhaven experiment E821: Experimental Technique 

In general terms, such as the use of electrostatic focusing and the magic 7, E821 
was modeled closely on the third CERN experiment. However, the statistical error 
goal and the ability of the Brookhaven AGS to provide the (necessary) tremendous 
increase in beam flux posed enormous challenges to the injection system as well as to 
the detectors, electronics and data acquisition. At the same time, systematic errors, 
in the measurements of both the field and anomalous precession frequencies, had to 
to be reduced significantly compared the the third CERN experiment. The challenges 
included: 

• How to allow multiple beam bunches from the AGS to pass undeflected through 
the back leg of the storage ring magnet. The use of multiple bunches reduces the 
instantaneous rates in the detectors, but it was judged technically impractical to 
rapidly pulse an infiector like that used at CERN. 

• How to provide a homogeneous storage ring field, uniform to 1 ppm when averaged 
over azimuth. 

• How to store an injected muon beam, which would greatly reduce the hadronic 
flash seen by the detectors. 
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• How to maximize the statistical power of the detected electrons 

• How to make good measurements on the decay electrons - good energy resolution, 
and, in the face of high data rates with a wide dynamic range, how to maintain the 
stability of signal gain and timing pickoff 

Some of the solutions were the natural product of technical advances over the 
previous twenty years. Others, such as the superconducting inflector, were truly novel, 
and required considerable development. A comparison of E821 and the third CERN 
experiment is given in Table 2. 

Table 2. A comparison of the features of the E821 and the third CERN muon 
{g — 2) experiment [38]. Both experiments operated at the "magic" 7 = 29.3, and 
used electrostatic quadrupolcs for vertical focusing. Bailey, et al. [38], do not quote a 
systematic error on the muon frequency Ua 



(Quantity 




CEKiN 


Magnet 


superconducting 


Room Temperature 


Yoke Construction 


Monolithic Yoke 


40 Separate Magnets 


Magnetic Field 


1.45 T 


1.47 T 


Magnet Gap 


180 mm 


140 mm 


Stored Energy 


6 MJ 




Field mapped in situ? 


yes 


no 


Central Orbit Radius 


7112 mm 


7000 mm 


Averaged Field Uniformity 


±1 ppm 


±10 ppm 


Muon Storage Region 


90 mm Diameter Circle 


120 X 80 mm^ Rectangle 


Injected Beam 


Muon 


Pion 


Inflector 


Static Superconducting 


Pulsed Coaxial Line 


Kicker 


Pulsed Magnetic 


IT ^ fi i/f^ decay 


Kicker Efficiency 


-4% 


125 ppm 


Muons stored/fill 


10^ 


350 


Ring Symmetry 


Four-fold 


Two-fold 


■\//3max/ Pm'm 


1.03 


1.15 


Detectors 


Pb-Scintillating Fiber 


Pb-Scintillator "Sandwich" 


Electronics 


Waveform Digitizers 


Discriminators 


Systematic Error on B-field 


0.17 ppm 


1.5 ppm 


Systematic Error on tOa 


0.21 ppm 


Not given 


Total Systematic Error 


0.28 ppm 


1.5 ppm 


Statistical Error on uja 


0.46 ppm 


7.0 ppm 


Final Total Error on 


0.54 ppm 


7.3 ppm 



An engineering run with pion injection occurred in 1997, a brief //-injection run 
where the new muon kicker was commissioned happened in 1998, and major data- 
collection periods of 3 to 4 months duration took place in 1999, 2000 and 2001. In 
each period, protons were accelerated by a hnear accelerator, accelerated further by a 
booster synchrotron, and then injected into the BNL alternating gradient synchrotron 
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(AGS). Radio- frequency cavities in the AGS ring provide acceleration to a momentum 
of 24 GeV/c, and maintain the protons in a number of discrete, equally spaced bunches. 
The number of bunches (harmonic number) in the AGS during a 2.7 s acceleration cycle 
was different for each of these periods: eight in 1999, six in 2000 and twelve in 2001. 
The AGS has the ability to dehver up to 70 x 10^^ protons (70 Tp) in one AGS cycle, 
providing a proton intensity per hour 180 times greater than that available at CERN in 
the 1970s. 



AGS 

U-Vline\ Xfc^^ 




Figure 4. The E821 bcamlinc and storage ring. Pions produced at 0° are collected 
by the quadrupoles Q1-Q2 and the momentum is selected by the collimators K1-K2. 
The pion decay channel is 72 m in length. Forward muons at the magic momentum 
are selected by the collimators K3-K4. 

The proton beam is extracted from the AGS one bunch at a time at 33 ms intervals. 
Each proton bunch results in a narrow time bunch of muons which is injected into 
the storage ring, and then the electrons from muon decays are measured for about 
10 muon lifetimes, or about 640 /is. A plan view of the Brookhaven Alternating 
Gradient Synchrotron, injection line and storage ring are shown in Figure 4. Because 
the maximum total intensity available from the AGS is < 70 Tp, the bunch intensity and 
the resulting pile-up (accidental coincidences between two electrons) in the detectors is 
minimized by maximizing the number of proton bunches. Pulse pile-up in the detectors 
following injection into the storage ring is one of the systematic issues requiring careful 
study in the data analysis. 
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2.1. The Proton and Muon Beamlines 

The primary proton beam from the AGS is brought to a water-cooled nickel production 
target. Because of mechanical shock considerations, the intensity of a bunch is limited 
to less than 7 Tp. The beamline, shown in Figure 4, accepts pions produced at 0° 
at the production target. They are collected by the first two quadrupoles, momentum 
analyzed, and brought into the decay channel by four dipoles. A pion momentum of 
3.115 GeV/c, 1.7 % higher than the magic momentum, is selected. The beam then enters 
a straight 80-meter-long focusing-defocusing quadrupole channel, where those muons 
from pion decays that are emitted approximately parallel to the pion momentum, so- 
called forward decays, are collected and transported downstream. Muons at the magic 
momentum of 3.094 MeV/c, and having an average polarization of 95%, are separated 
from the slightly higher momentum pions at the second momentum slit. However, 
after this momentum selection a rather large pion component remains in the beam, 
whose composition was measured to be 1:1:1, e"*" : /i"*" : tt"^. The proton content was 
calculated to be approximately one-third of the pion flux [48]. The secondary muon 
beam intensity incident on the storage ring was about 2 x 10^ per fill of the ring, which 
can be compared with 10® particles per fill with "pion injection[47]" which was used in 
the 1997 engineering run. 

The detectors are exposed to a large flux of background pions, protons, neutrons, 
electrons and other particles associated with the beam injection process, called the 
injection "flash". The intensity varies around the ring, being most intense in the 
detectors adjacent to the injection point (referred to below as 'upstream' detectors). The 
flash consists of prompt and delayed components, with the prompt component caused by 
injected particles passing directly into the detectors. The delayed component is mainly 
due to 7 rays following neutron capture. The neutrons are produced primarily by the 
protons and pions in the beam striking the magnet, calorimeters, etc. Many neutrons 
thermalize in the calorimeter and surrounding materials over tens of microseconds, where 
they can undergo nuclear capture. These 7 rays cause a DC baseline offset in the PMT 
signal which steadily decays with a lifetime of ~ 50-100 yus. To reduce the decay time 
of the neutron background, the epoxy in the upstream subset of detectors is doped 
with natural boron carbide powder, thus taking advantage of the large neutron capture 
cross-section on ^^B. 

This injection flash is most severe with the pion injection scheme. The upstream 
photomultiplier tubes had to be gated off for 120 /is (1.8 muon lifetimes) following 
injection, to allow the signals to return to the nominal baseline. To reduce this "flash" 
and to increase significantly the number of muons stored per fill of the storage ring, a 
fast muon kicker was developed which permitted direct muon injection into the storage 
ring. 
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2.2. The Inflector Magnet 

The geometry of the incoming beam is shown in Figure 5. A unique superconducting 
inflector magnet [58, 59] was built to cancel the magnetic field permitting the beam to 
arrive undeflected at the edge of the storage region (see Figure 6(b)). The inflector is 
a truncated, double-cosine theta magnet, shown in Figure 5(b) at its downstream end, 
with the muon velocity going into the page. In the inflector, the current flows into the 
page down the central "C" -shaped layer of superconductor, then out of the page through 
the "backward-D" -shaped outer conductor layer. At the inflector exit, the center of the 
injected beam is 77 mm from the central orbit. For /i"*" stored in the ring, the main field 
points up, and qv x B points to the right in Figure 5(b), toward the ring center. 



R = 71 12 mm from ring center 



- Tangential reference line 





Beam vacuum 
clramber 



Superconducting \ -Partition wall 
coils Passive superconducting 

shield 



(a) 



(b) 



Figure 5. (a) A plan view of the inflector-storage-ring geometry. The dot-dash 
hne shows the central muon orbit at 7112 mm. The beam enters through a hole in 
the back of the magnet yoke, then passes into the inflector. The inflector cryostat 
has a separate vacuum from the beam chamber, as can be seen in the cross sectional 
view. The cryogenic services for the inflector are provided through a radial penetration 
through the yoke at the upstream end of the inflector. (b) A cross-sectional view of 
the pole pieces, the outer-radius coil-cryostat arrangement, and the downstream end 
of the superconducting inflector. The muon beam direction at the inflector exit is into 
the page. The center of the storage ring is to the right. The outer-radius coils which 
excite the storage-ring magnetic field are shown, but the inner-radius coils are omitted. 



In this design shown in Figure 5(b) and Figure 6(b), the beam channel aperture 
is rather small compared to the flux return area. The field is ~ 3 T in the return 
area (inflector plus central storage- ring field), and the fiux density is sufficiently high to 
lower the critical current in superconductor placed in that region. If the beam channel 
aperture were to be increased by pushing the coil further into the flux return area, 
the design would have to be changed, either by employing a superconductor with larger 
critical current, or by using more conductors in a revised geometry, further complicating 
the fabrication of this magnet. The result of the small inflector aperture is a rather poor 
phase-space match between the inflector and the storage ring and, as a consequence, a 
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loss of stored muons. 

As can be seen from Figure 6(a), the entrance (and exit) to the beam channel is 
covered with superconductor, as well as by aluminum windows that are not visible in 
the photograph. This design was chosen to maximize the mechanical stability of the 
superconductor in the magnetic field, thus reducing the risk of motion which would 
quench the magnet. However, multiple scattering in the material at both the entrance 
and exit windows causes about half the incident muon beam to be lost. 

The distribution of conductor on the outer surface of the inflector magnet (the 
"D-shaped" arrangement) prevents most of the magnetic flux from leaking outside of 
the inflector volume, as seen from Figure6(b). To prevent flux leakage from entering 
the beam storage region, the inflector is wrapped with a passive superconducting shield 
that extends beyond both inflector ends, with a 2 m seam running longitudinally along 
the inflector side away from the storage region. With the inflector at zero current and 
the shield warm, the main storage ring magnet is energized. Next the inflector is cooled 
down, so that the shield goes superconducting and pins the precision field inside the 
inflector region. When the inflector magnet is powered, the supercurrents in the shield 
prevent the leakage field from penetrating into the storage region behind the shield. 




(a) (b) 

Figure 6. (a) A photo of the prototype inflector showing the crossover between the 
two coils. The beam channel is covered by the lower crossover, (b) The magnetic 
design of the inflector. Note that the magnetic flux is largely contained inside of the 
inflector volume. 

A 0.5 m long prototype of the inflector was fabricated and then tested, first 
by itself, then in an external magnetic field. Afterward the full-size 1.7 m inflector 
was fabricated. Unfortunately, the full-size inflector was damaged during initial tests, 
making it necessary to cut through the shield to repair the interconnect between the 
two inflector coils. A patch was placed over the cut, and this repaired inflector was used 
in the 1997 7r-injection run as well as in the 1998 and 1999 run periods. Magnetic flux 
leakage around the patch reduced the main storage-ring field locally by ~ 600 ppm. 
This inhomogeneity complicated the field measurement due to the degradation in NMR 
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performance, and contributed an additional uncertainty of 0.2 ppm to the average 
magnetic field seen by the muons. The leakage field from a new inflector, installed 
before the 2000 running period, was immeasurably small. 

2.3. The Fast Muon Kicker 

Left undisturbed, the injected muon beam would pass once around the storage ring, 
strike the inflector and be lost. As shown in Figure 7, the role of the fast muon kicker 
is to briefly reduce a section of the main storage field and move the center of the muon 
orbit to the geometric center of the storage ring. The 77-mm offset at the injection 
point, between the center of the entering beam and the central orbit, requires that 
the beam be kicked outward by approximately 10 mrad. The kick should be made at 
about 90 degrees around the ring, plus a few degree correction due to the defocusing 
effect of the electric quadrupoles between the injection point and the kicker, as shown 
in Figures 5(b) and 7. 



Figure 7. A sketch of the beam geometry. R — 7112 mm is the storage ring radius, 
Xc = 77 mm is the distance between the inflector center and the center of the storage 
region. This is also the distance between the centers of the circular trajectory that 
a particle entering at the inflector center (at = with x' = 0) will follow, and the 
circular trajectory a particle at the center of the storage volume (at 6* = with x' = 0) 
will follow. 



The requirements on the fast muon kicker are rather stringent. While electric, 
magnetic, and combination electromagnetic kickers were considered, the collaboration 
settled on a magnetic kicker design [72] because it was thought to be technically easier 
and more robust than the other options. Because of the very stringent requirements 
on the storage ring magnetic field uniformity, no magnetic materials could be used. 
Thus the kicker field had to be generated and shaped solely with currents, rather than 
using ferrite cores. Even with the kicker field generated by currents, there existed the 
potential problem of inducing eddy currents which might affect the magnetic field seen 
by the stored muons. 

The length of the kicker is limited to the ~ 5 m azimuthal space between the 
electrostatic quadrupoles (see Figure 10), so each of the three sections is 1.76 m long. 
The cross section of the kicker is shown in Figure 8. The two parallel conductors are 
connected with cross-overs at each end, forming a single current loop. The kicker plates 




n 
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also have to serve as "rails" for the NMR field- mapping trolley (discussed below), and 
the trolley is shown riding on the kicker rails in Figure 8. 

Macor HV Standoff 




Macor Baseplate 

Figure 8. An elevation view of the kicker plates, showing the ceramic cage supporting 
the kicker plates, and the NMR trolley riding on the kicker plates. (The trolley is 
removed during data collection). 

The kicker current pulse is formed by an under-damped LCR circuit. A capacitor is 
charged to 95 kV through a resonant charging circuit. Just before the beam enters the 
storage ring, the capacitor is shorted to ground by firing a deuterium thyratron. The 
peak current in an LCR circuit is given by /q = VQ/{uJdL) making it necessary to keep 
the system inductance, L, low to maximize the magnetic field for a given voltage Vq. 
For this reason, the kicker was divided into three sections, each powered by a separate 
pulse-forming network. 

The electrode design and the time-varying fields were modeled using the package 
OPERA 2d [73]. Calculations for electrodes made of Al, Ti and Cu, were carried out with 
the geometry shown in Figure 8. The aluminum electrodes were best for minimizing eddy 
currents 20 /is after injection, which was fortunate, since mechanical tests on prototype 
electrodes made from the three metals showed that only aluminum was satisfactory. 
This was also true from the detector point of view, since the low atomic number of 
aluminum minimizes the multiple scattering and showering of the decay electrons on 
their way to the calorimeters. The magnetic field produced by the kicker, along with 
the residual field, was measured using the Faraday effect [72], and these are shown in 
Figure 9. 

The cyclotron period of the ring is 149 ns, substantially less than the kicker pulse 
base-width of ~ 400 ns, so that the injected beam is kicked on the first few turns. 
Nevertheless, approximately 10^ muons are stored per fill of the ring, corresponding to 
an injection efficiency of about 3 to 5% (ratio of stored to incident muons). The storage 
efficiency with muon injection is much greater than that obtained with pion injection 
for a given proton fiux, with only 1 percent of the hadronic fiash. 
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Figure 9. (a)Thc main magnetic field of tlie kicker measured -with the Faraday effect, 
(b) The residual magnetic field measured using the Faraday effect. The solid points are 
calculations from OPERA, and the small x are the experimental points measured ■with 
the Faraday effect. The solid horizontal lines sho^w the ±0.1 ppm band for affecting 
J B-dl 

2.4- The Electrostatic Quadrupoles 

The electric quadrupoles, •which are arranged around the ring •with four-fold symmetry, 
provide vertical focusing for the stored muon beam. The quadrupoles cover 43% of the 
ring in azimuth, as shown in Figure 10. While the ideal vertical profile for a quadrupole 
electrode would be hyperbolic, beam dynamics calculations determined that the higher 
multipoles present with flat electrodes, which are much easier to fabricate, would not 
cause an unacceptable level of beam losses. The flat electrodes are shown in Figure 11. 
Only certain multipoles are permitted by the four-fold symmetry, and a judicious choice 
of the electrode width relative to the separation between opposite plates minimizes 
the lowest of these. With this configuration, the 20-pole is the largest, being 2% of 
the quadrupole component and an order of magnitude greater than the other allowed 
multipoles [74]. 

In the quadrupole regions, the combined electric and magnetic fields can lead 
to electron trapping. The electron orbits run longitudinally along the inside of the 
electrode, and then return on the outside. Excessive trapping in the relatively modest 
vacuum of the storage ring can cause sparking. To minimize trapping, the leads were 
arranged to introduce a dipole field at the end of the quadrupole thus sweeping away 
trapped electrons. In addition, the quadrupoles are pulsed, so that after each fill of the 
ring all trapped particles could escape. Since some protons (antiprotons) were stored in 
each fill of /i"*" they were also released at the end of each storage time. This lead 

arrangement worked so well in removing trapped electrons that for the /i"*" polarity it 
would have been possible to operate the quadrupoles in a dc mode. For the storage of 
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Inflector 




Figure 10. The layout of the storage rmg, as seen from above, showmg the location 
of the inflector, the kicker sections (labeled K1-K3), and the quadrupoles (labeled Ql- 
Q4). The beam circulates in a clockwise direction. Also shown are the collimators, 
which are labeled "C", or "^0" indicating whether the Cu collimator covers the full 
aperture, or half the aperture. The collimators are rings with inner radius: 45 mm, 
outer radius: 55 mm, thickness: 3 mm. The scalloped vacuum chamber consists of 12 
sections joined by bellows. The chambers containing the inflector, the NMR trolley 
garage, and the trolley drive mechanism are special chambers. The other chambers 
are standard, with either quadrupolc or kicker assemblies installed inside. An electron 
calorimeter is placed behind each of the radial windows, at the postion indicated by 
the calorimeter number. 

this was not true; the trapped electrons necessitated an order of magnitude better 
vacuum, and hmited the storage time to less than 700 yus. 

Beam losses during the measurement period, which could distort the expected time 
spectrum of decay electrons, had to be minimized. Beam scraping is used to remove, 
just after injection, those muons which would likely be lost later on. To this end, 
the quadrupoles are initially powered asymmetrically, and then brought to their final 
symmetric voltage configuration. The asymmetric voltages lower the beam and move it 
sideways in the storage ring. Particles whose trajectories reach too near the boundaries 
of the storage volume (defined by collimators placed at the ends of the quadrupole 
sectors) are lost. The scraping time was 17 /is during all data collection runs except 
2001, where 7 yus was used. The muon loss rates without scraping were on the order of 
0.6% per lifetime at late times in a fill, which dropped to ~ 0.2% with scraping. 

2.5. The Superconducting Storage Ring 

The storage ring magnet combined with the electrostatic quadrupoles form a Penning 
trap that, while very different in scale, has common features with the electron ^f-value 




Figure 11. A photograph of an electrostatic quadrupole assembly inside a vacuum 
chamber. The cage assembly doubles as a rail system for the NMR trolley which is 
resting on the rails. The location of the NMR probes inside the trolley are shown 
as black circles (see §2.6). The probes arc located just behind the front face. The 
inner (outer) circle of probes has a diameter of 3.5 cm (7 cm) at the probe centers. 
The storage region has a diameter of 9 cm. The vertical location of three upper fixed 
probes is also shown. The fixed probes are located symmetrically above and below the 
vacuum chamber. 

measurements [21, 54]. In the electron experiments, a single electron is stored for long 
periods of time, and the spin and cyclotron frequencies are measured. Since the muon 
decays in 2.2 yus , the study of a single muon is impossible. Rather, an ensemble of muons 
is stored at relativistic energies, and and their spin precession in the magnetic field is 
measured using the parity- violating decay to analyze the spin motion(see Equations 17, 
18, 8). 

The design goal of ±1 ppm was placed on the field uniformity when averaged over 
azimuth in the storage ring. A "superferric" design, where the field configuration is 
largely determined by the shape and magnetic properties of the iron, rather than by the 
current distribution in the superconducting coils, was chosen. To reach the ppm level of 
uniformity it was important to minimize discontinuities such as holes in the yoke, spaces 
between adjacent pole pieces, and especially the spacing between pole pieces across the 
magnet gap containing the beam vacuum chamber. Every effort was made to minimize 
penetrations in the yoke, and where they are necessary, such as for the beam entrance 
channel, additional iron is placed around the hole to minimize the effect of the hole on 
the magnetic flux circuit. 

The storage ring, shown in Figure 12, is designed as a continuous C-magnet[60] 
with the yoke made up of twelve sectors with minimum gaps where the yoke pieces 
come together. A cross-section of the magnet is shown in Figure 13. The largest gap 
between adjacent yoke pieces after assembly is 0.5 mm. The pole pieces are built in 
36 pieces, with keystone rather than radial boundaries to ensure a close fit. They 
are electrically isolated from each other with 80 /im kapton to prevent eddy currents 
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from running around the ring, especially during a quench or energy extraction from the 
magnet. The vertical mismatch from one pole piece to the next when going around the 
ring in azimuth is held to ±10 /im, since the field strength depends critically on the 
pole-piece spacing across the magnet gap. 




Figure 12. The storage-ring magnet. Tiie cryostats for tlie inner-radius coils are 
clearly visible. The kickers have not yet been installed. The racks in the center are 
the quadrupole pulsers, and a few of the detector stations are installed, especially the 
quadrant of the ring closest to the person. The magnet power supply is in the upper 
left, above the plane of the ring. (Courtesy of Brookhaven National Laboratory) 

The field is excited by 14 m-diameter superconducting coils, which in 1996 were 
the largest-diameter such coils ever fabricated. The coil at the outer radius consists of 
two identical coils on a common mandrel, above and below the plane of the beam, each 
with 24 turns. Each of the inner-radius coils, which are housed in separate cryostats, 
also consists of 24 turns (see Figures 5(b) and 13(a)). The nominal operating current 
is 5200 A, which is driven by a power supply. The choice of using an extremely stable 
power supply, further stabilized with feedback from the NMR system, was chosen over 
operating in a "persistent mode" , for two reasons. The switch required to change from 
the powering mode to persistent mode was technically very complicated, and unlike the 
usual superconducting magnet operated in persistent mode, we anticipated the need to 
cycle the magnet power a number of times during a three-month running period. 

The pole pieces are fabricated from continuous vacuum-cast low-carbon magnet 
steel (0.0004% carbon), and the yoke from standard AISI 1006 (0.07% carbon) magnet 
steel[60]. At the design stage, calculations suggested that the field could be made 
quite uniform, and that when averaged over azimuth, a uniformity of ±1 ppm could be 
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Figure 13. (a) A Cross-sectional view of the storage ring magnet. The beam center 
is at a radius of 7112 mm. The pole pieces are separated from the yoke by an air gap. 
(b) An expanded view of the pole pieces, which shows the edge shims and the wedge 
shim. 

achieved. It was anticipated that, at the initial turn-on of the magnet, the field would 
have a uniformity of about 1 part in 10^, and that an extensive program of shimming 
would be necessary to reach a uniformity of one ppm. 

A number of tools for shimming the magnet were therefore built into the design. 
The air gap between the yoke and pole pieces dominates the reluctance of the magnetic 
circuit outside of the gap that includes the storage region, and decouples the field in 
the storage region from possible voids, or other defects in the yoke steel. Iron wedges 
placed in the air gap were ground to the wedge angle needed to cancel the quadrupole 
field component inherent in a C- magnet. The dipole can be tuned locally by moving the 
wedge radially. The edge shims screwed to the pole pieces were made oversize. Once the 
initial field was mapped and the sextupole content was measured, the edge shims were 
reduced in size individually in two steps of grinding, followed by field measurements, to 
cancel the sextupole component of the field. 

After mechanical shimming, these higher multipoles were found to be quite constant 
in azimuth. They are shimmed out on average by adjusting currents in conductors 
placed on printed circuit boards going around the ring in concentric circles spaced by 
0.25 cm. These boards are glued to the top and bottom pole faces between the edge shims 
and connected at the pole ends to form a total of 240 concentric circles of conductor, 
connected in groups of four, to sixty ±1 A power supplies. These correction coils are 
quite effective in shimming multipoles up through the octupole. Multipoles higher than 
octupole are less than 1 ppm at the edge of the storage aperture, and, with our use of 
a circular storage aperture, are unimportant in determining the average magnetic field 
seen by the muon beam (see § 4.1.4). 
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2.6. Measurement of the Precision Magnetic Field 

The magnetic field is measured and monitored by pulsed Nuclear Magnetic Resonance 
tecliniques[75]. The free-induction decay is picked up by the coil Ls in Figure 14 after 
a pulsed excitation rotates the proton spin in the sample by 90° to the magnetic field, 
and is mixed with the reference frequency to form the frequency /fid- The reference 
frequency, /ref = 61.74 MHz, is obtained from a frequency synthesizer phase locked 
to the LORAN C standard[77], that is chosen with /rcf < /nmr such that for all 
probes /fid — 50 kHz. The relationship between the actual field i?reai and the field 
corresponding to the reference frequency is given by 

-Breal = -Bfcf ( 1 + ^^^] ■ (22) 
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Figure 14. Schematics of different NMR probes, (a) Absolute probe featuring a 
spherical sample of water. This probe was the same one used in reference [39] to 
determine A, the muon-to-proton magnetic-moment ratio, (b) Plunging probe, which 
can be inserted into the vacuum at a specially shimmed region of the storage ring 
to transfer the calibration to the trolley probes, (c) The standard probes used in 
the trolley and as fixed probes. The resonant circuit is formed by the two coils with 
inductances Lg and Lp and a capacitance Cg made by the Al-housing and a metal 
electrode. 

The field measurement process has three aspects: calibration, monitoring the field 
during data collection, and mapping the field. The probes used for these purposes are 
shown in Figure 14. To map the field, an NMR trolley was built with an array of 17 
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NMR probes arranged in concentric circles, as shown in Figure 11. While it would be 
preferable to have information over the full 90-mm aperture, space limitations inside 
the vacuum chamber, which can be understood by examining Figures 8 and 11, prevent 
a larger diameter trolley. The trolley is pulled around the storage ring by two cables, 
one in each direction circling the ring. One of these cables is a standard thin (Lemo- 
connector compatible) co-axial signal cable, and the other is non-conducting because 
of its proximity to the kicker high voltage. Power, data in and out, and a reference 
frequency all are carried on the single Lemo cable. 

During data-collection periods, the trolley is parked in a garage (see Figure 10) 
in a special vacuum chamber. Every few days, at random times, the field is mapped 
using the trolley. To map the field, the trolley is moved into the storage region and 
pulled through the vacuum chamber, measuring the field at over 6000 locations. During 
the data collection periods, the storage-ring magnet remains powered continuously for 
periods lasting from five to twenty days, thus the conditions during mapping are identical 
to those during the data collection. 

To cross calibrate the trolley probes, bellows are placed at one location in the 
ring to permit a special NMR plunging probe, or a absolute calibration probe with a 
spherical water sample[76], to plunge into the vacuum chamber. Measurements of the 
field at the same spatial point with the plunging, calibration and trolley probes provides 
both relative and absolute calibration of the trolley probes. During the calibration 
measurements before and after each running period, the spherical water probe is used 
to calibrate the plunging probe, the center probe of the trolley, and several other trolley 
probes. The absolute calibration probe provides the calibration to the Larmor frequency 
of the free proton[65], which is called Up below. 

To monitor the field on a continuous basis during data collection, a total of 378 
NMR probes are placed in fixed locations above and below the vacuum chamber around 
the ring. Of these, about half provide data useful in monitoring the field with time. 
Some of the others are noisy, or have bad cables or other problems, but a significant 
number of fixed probes are located in regions near the pole-piece boundaries where 
the magnetic gradients are sufficiently large to reduce the free-induction decay time 
in the probe, limiting the precision on the frequency measurement. The number of 
probes at each azimuthal position around the ring alternates between two and three, at 
radial positions arranged symmetrically about the magic radius of 7112 mm. Because 
of this geometry, the fixed probes provide a good monitor of changes in the dipole and 
quadrupole components of the field around the storage ring. 

Initially the trolley and fixed probes contained a water sample. Over the course of 
the experiment, the water samples in many of the probes were replaced with petroleum 
jelly. The jelly has several advantages over water: low evaporation, favorable relaxation 
times at room temperature, a proton NMR signal almost comparable to that from water, 
and a chemical shift (and the accompanying NMR frequency shift) with a temperature 
coefficient much smaller than that of water, and thus negligible for our experiment. 
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2. 7. Detectors and Electronics 

The detector system consists of a variety of particle detectors: calorimeters, position- 
sensitive hodoscope detectors, and a set of tracking chambers. There are also horizontal 
and vertical arrays of scintillating-fiber hodoscopes which could be temporarily 
inserted into the storage region. A number of custom electronics modules were 
developed, including event simulators, multi-hit time-to-digital converters(MTDC), and 
the waveform digitizers (WFD) which are at the heart of the measurement. We refer 
to the data collected from one muon injection pulse in one detector as a 'spill,' and we 
will speak of 'early-to-late' effects: namely, the gain or time stability requirements in a 
given detector at early compared to late decay times in a spill. 

2.7.1. Electromagnetic Calorimeters: Design Considerations The electromagnetic 
calorimeters, together with the custom WFD readout system, are the primary source 
of data for determining the precession frequency. They provide the energies and arrival 
times of the electrons, and they also provide signal information immediately before and 
after the electron pulses, allowing studies of baseline changes and pulse pile-up. 

There are 24 calorimeters placed evenly around the 45-m circumference of the 
storage region, adjacent to the inside radius of the storage vacuum region as shown 
in Figures 10 and 12. Nearly all decay electrons have momenta (0 < p{lab) < 3.1 
GeV/c, see Figure 2) below the stored muon momentum (3.1 GeV/c ±0.2%), and they 
are swept by the B-field to the inside of the ring where they can be intercepted by the 
calorimeters. The storage-region vacuum chamber is scalloped so that electrons pass 
nearly perpendicular to the vacuum wall before entering the calorimeters, minimizing 
electron pre-showering (see Figure 10). The calorimeters are positioned and sized in 
order to maximize the acceptance of the highest-energy electrons, which have the largest 
statistical figure of merit NA"^. The variations of and A as a function of electron 
energy are shown in Figures 2 and 3. 

The electrons with the lowest laboratory energies, while more numerous than high 
energy electrons, generally have a lower figure of merit and therefore carry relatively 
little information on the precession frequency. These electrons have relatively small 
radii of curvature, and exit the the ring vacuum chamber closer to the radial direction 
than electrons at higher energies, with most of them missing the detectors entirely. 
Detection of these electrons would require detectors that cover a much larger portion of 
the circumference than is needed for high-energy electrons, and is not cost effective. 

Consequently the detector system is designed to maximize the acceptance of the 
high energy decay electrons above approximately 1.8 GeV, with the acceptance falling 
rapidly below this energy. The detector acceptance reaches a maximum of 87% at 2.3 
GeV, decreases to 70% at 1.8 GeV, and continues to decrease roughly linearly to zero 
as the energy decreases. With increasing energy above 2.3 GeV, the acceptance also 
decreases because the highest-energy electrons tend to enter the calorimeters at the 
outer radial edge, increasing the loss of registered energy due to shower leakage, and 
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reducing the acceptance to 80% at 3 GeV. 

In a typical analysis, the full data sample consists of all electrons above a threshold 
energy of about 1.8 GeV, where NA^ is approximately a maximum, with about 65% 
of the electrons above that energy detected (Figure 3). The average asymmetry 
is about 0.35. The loss of efficiency is from the low-energy tail in the detector 
response characteristic of electromagnetic showers in calorimeters, and from lower energy 
electrons missing the detectors altogether. The statistical error improves by only 5% if 
the data sample contains all electrons above 1.8 GeV compared to all above 2.0 GeV. 
For threshold energies below 1.7 GeV, the decline in the average asymmetry more than 
cancels the additional number of electrons in NA^, and the statistical error actually 
increases. Some of the independent analyses fit time spectra of data formed from 
electrons in narrow energy bands (about 200 MeV wide). When the results of the 
separate fits are combined, there was is a 10% reduction in the statistical error on uja- 
however, there is also a slight increase in the systematic error contribution from gain 
shifts, because the relative number of events moved by a gain shift from one energy band 
to another increased. One analysis used data weighted by the asymmetry as a function 
of energy. It can be shown[71] that this produces the same statistical improvement as 
dividing the data into energy bands. 

Gain and timing shift limitations are much more stringent within a single spill than 
from spill to spill. Shifts at late decay times compared to early times in a given spill, 
so-called 'early-to-late' shifts, can lead directly to serious systematic errors on uja- Shifts 
of gain or the t = point from one spill to the next are generally much less serious; they 
will usually only change the asymmetry, average energy, phase, etc., hut to first order, 
ijja will be unaffected. 

The calorimeters should have pulses with narrow time widths to minimize the 
probability of two pulses overlapping (pile-up) during the very high electron decay data 
rates encountered at early decay times, which can reach a MHz in a single detector. The 
scintillator is chosen to have minimal long-lived components to reduce the afterglow from 
the intense detector fiash associated with beam injection. Laser calibration studies show 
that the timing stability for a typical detector over any 200 microsecond time interval is 
better than 15 ps, easily meeting the demands of the measurement of a^. For example, 
a 20 ps timing shift would lead to an uncertainty in of about 0.1 ppm, which is small 
compared to the final error. Modest detector energy resolution (^ 10 — 15% at 2 GeV) is 
required in order to select the desired high energy electrons for analysis. Better energy 
resolution also reduces the amount of calibration data needed to monitor the stability 
of the detector gains. 

The stability requirement for the electron energy measurement ('gain') versus time 
in the spill is largely determined by the energy dependence of the phase of the {g — 2) 
oscillation. In a fit of the data to the 5-parameter function, the oscillation phase is 
highly correlated to Ua- Therefore a shift in the gain from early to late decay times, 
combined with an energy dependence in the (g-2) phase, can lead to a systematic error 
in the determination of uJa- There are two main contributing factors to the energy 
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dependence, which appear with opposite signs: 1) The phase in the 5-parameter 
function (Equation 8) depends on the electron drift time. High-energy electrons must 
travel further, on average, from the point of muon decay to the detector and therefore 
have longer drift times. The change in drift time with energy implies a corresponding 
energy dependence in the {g — 2) phase. 2) For decay electrons at a given energy, 
those with positive (radially outward) components of momentum at the muon decay 
point travel further to reach the detectors than electrons with negative (radially inward) 
components. They spread out more in the vertical direction and may miss the detectors 
entirely. Consequently, electrons with positive (outward) radial momentum components 
will have slightly lower acceptance than those with negative components, causing the 
average spin direction to rotate slightly, leading to a shift in 0. Recalling the correlation 
between electron direction and muon spin, the overall effect is to shift the time when 
the number oscillation reaches its maximum, causing a shift in the precession phase in 
the 5-parameter function. The size of the shift depends on the electron energy. From 
studies of the data sample and simulations, it is established that the detector gains need 
to be stable to better than 0.2% over any 200 /is time interval in a spill, in order to 
keep the systematic error contribution to uja less than 0.1 ppm from gain shifts. This 
requirement is met by all of the calorimeters. The gain from one spill to the next is 
not coupled to the the precession frequency and therefore the requirement on the spill 
to spill stability is far less stringent than the stability requirement within an individual 
spill. 

2.7.2. Electromagnetic Calorimeters: Construction The calorimeters (see Figure 15) 
consist of 1 mm diameter scintillating fibers embedded in a lead-epoxy matrix. [61] 
The final dimensions are 22.5 cm (radial) x 14 cm high x 15 cm deep, where the 
vertical dimension is limited by the vertical magnet gap. The radiation length is 
Xq = 1.14 cm with a Moliere radius of about 2.5 cm. This gives a calorimeter depth 
of 13 radiation lengths, which means between 96% and 93% of the energies from 1 
to 3 GeV are deposited in the calorimeter, with the remaining energy escaping. The 
fibers are oriented along the radial direction in the ring, so that incoming electrons 
are incident approximately perpendicular to the fibers. The large-radius ends of the 
fibers are reflective in order to improve the uniformity of the pulse height response as a 
function of the radial direction. There is an average 6% change in the photomultiplier 
(PMT) pulse height depending on the radial electron entrance point, and there is a 
change of a few percent depending on the vertical position. The measured versus actual 
electron energy deviates from a linear relationship by about 1% between and 3 GeV, 
primarily due to shower leakage. 

Of the electron shower energy deposited in the calorimeter, about 5% is deposited 
in the scintillator, producing on the order of 500 photo-electrons per GeV of incident 
electron energy, with the remainder being deposited in the lead or epoxy. The average 
overall energy resolution is a{E)/E ^ 10%/ ^E{GeV), dominated by the statistics of 
shower sampling (i.e. variation in the energy that gets deposited in the scintillator vs. 
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Figure 15. Schematic of a detector station. The calorimeter consists of scintillating 
fibers embedded in lead, with the fibers being directed radially toward the center of 
the storage ring. The electron entrance face of the calorimeter is covered with an array 
of five horizontally-oriented FSD scintillators. Six of the detector stations also have a 
PSD scintillating xy hodoscope covering the entrance face. 

the other calorimeter material), with some contribution from photo-electron statistics 
and shower leakage out the back or sides of the calorimeter. 

Scintillation light is gathered by four 52.5 cm^, nearly square, acrylic light-guides 
which cover the inner radius side of the calorimeter. Each light-guide necks down to a 
circle which mates to a 12-dynode 5 cm-diameter Hamamatsu R1828 photo- multiplier 
tube. The photomultiplier gains in the four segments of each calorimeter are balanced 
initially using monochromatic 2 GeV electrons from a BNL test beam. Some detector 
calibrations shifted when the detectors were moved to their final storage ring locations. 
In these cases the gains in the quadrants were balanced using the energy spectra from 
the muon decay electrons themselves. The electron signals, which are the sum of the 4 
signals from the photomulutipliers on the calorimeter quadrants, are about 5 ns wide 
(FWHM), and are sent to waveform digitizers for readout. 

In one spill during the time interval from a few tens of microseconds to 640 /is after 
injection, approximately 20 decay-electrons above 1.8 GeV are recorded on average 
in each detector. The instantaneous rate of decay electrons above 1 GeV changes from 
about 300 kHz to almost zero over this period. The gain and timing of photo-multipliers 
can depend on the data rate. The necessary gain and timing stability is achieved with 
custom, actively stabihzed photo-multiplier bases[62]. 

To prevent paralysis of the photo-multipliers due to the injection flash, the 
amplifications in the photo-multipliers are temporarily reduced by a factor of about 
1 million during the beam injection. Depending on the intensity of the flash and the 
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duration of the background levels encountered at a particular detector station, the 
amplifications are restored at times between 2 to 50 /is after injection. The switching 
of the gain is accomplished in the Hamamatsu R1828 photo-multipliers by swapping 
the bias voltages on dynodes 4 and 7. With the proper selection of the delay time 
after injection to let backgrounds die down, the gains typically return to 99.8% of their 
steady-state value within several microseconds after the tube is turned back on. Other 
gating schemes, such as switching the photocathode voltage, were found to have either 
required a much longer time for the gain to recover, or failed to give the necessary 
reduction in the gain when the tube is gated off. 

2.7.3. Special Detectors: FSD, PSD, Fiber Harps, Traceback Chambers Several 
specialized detector systems, the Front Scintillation Detectors (FSD), Position-sensitive 
Detectors (PSD), Fiber Harp Scintillators, and Traceback Detectors, are all designed 
to provide information on the phase-space parameters of the stored muon beam and 
their decay electrons. Such measurements are compared to simulation results, and are 
important, for example, in the study of coherent betatron motion of the stored beam 
and detector acceptances, and in placing a limit on the electric dipole moment of the 
muon (see § 3 and § 5). A modest knowledge of the beam phase space is necessary in 
order to calculate the average magnetic and electric fields seen by the stored muons. 

The FSDs cover the front (electron entrance) faces of the calorimeters at about 
half the detector stations. They provide information on the vertical positions of the 
electrons when they enter the calorimeters. Each consists of an array of five horizontally- 
oriented strips of scintillator 23.5 cm long (radial dimension) x 2.8 cm high x 1.0 cm 
thick. Each strip is coupled to a 2.8-cm-diameter Hamamatsu R6427 PMT. The 
signal is discriminated and sent to multi-hit time-to-digital converters (MTDC) for time 
registration. 

The FSDs are an integral part of the muon loss monitoring system. A lost muon 
is identified as a particle which passes through three successive detector stations, firing 
the three successive FSD detectors in coincidence. The candidate lost muon is also 
required to leave minimal energy in each of the calorimeters, since a muon loses energy 
by ionization, rather than creating an electromagnetic shower. This technique only 
identifies muons lost in the inward radial direction. Those lost in other directions, for 
example, above or below the storage region, are not monitored. The scaling between 
registered and actual muon losses is estimated based on models of the muon loss 
mechanisms. 

Five detector stations are equipped with PSDs. They provide better vertical 
position resolution than the FSDs, and also provide radial position information. Each 
PSD consists of an xy array of scintillator sticks which, like the FSDs, cover the electron 
entrance face of the calorimeters. There are 20 horizontally- and 32 vertically-oriented 
7 mm-wide, 8 mm-thick elements read out with multi-anode photomultipliers followed 
by discriminators and MTDCs. 

The traceback-chamber system consists of four sets of drift chambers which can 
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provide up to 12 vertical and 12 horizontal position measurements along a track, with 
about 300-micron resolution. It is positioned 270 degrees around the ring from the 
injection point, where a thin window is installed for the electrons to pass through the 
vacuum wall with minimal scattering. A regular calorimeter is placed behind the array 
to serve as a trigger. Using the momentum derived from the track information, plus 
a knowledge of the B-field, it is possible to extrapolate the track back to the point 
where the trajectory is tangent to the storage ring. Since the decay electrons are nearly 
parallel to the muon momentum, and the muons are traveling with small angles relative 
to the tangent, the extrapolation point is a very good approximation {ay ~ 9 mm, 
(7i? ~ 15 mm) to the position of muon decay. Thus, the traceback provides information 
on the muon distribution for that portion of the ring just upstream of the traceback 
position. 

The storage ring is equipped with two sets of scintillating-fiber beam monitors 
inside the vacuum chamber. Their purpose is to measure the beam profile as a function 
of time, and they can be inserted at will into the storage region. Each set consists of 
an X and y plane of seven 0.5-mm diameter scintillating fibers with 13 mm spacing that 
covers the central part of the beam region as shown in Figure 16. These fiber "harps" are 
located at 180° and 270° around the ring from the injection point as shown in Figure 10. 
Each fiber is mated to a clear fiber, which connects to a vacuum optical feed-through 
and then to a PMT which is read out with a WFD that digitizes continuously during the 
measurement time. The scintillating fibers are sufficiently thin that the beam can be 
monitored for many tens of microseconds after injection, as shown in Figure 21, before 
the stored beam is degraded (effective lifetime about 30 microseconds ), and therefore 
they are very helpful in providing information on the early-time position and width 
distributions of the stored muon beam. 




X monitor y monitor 

Figure 16. A sketch of the fiber beam monitors. Tlie fibers can be rotated into the 
beam for caUbration such that each fiber sees the same portion of the beam. 

2.7.4- Electronics For each calorimeter, the arrival times of the signals from the four 
photo-multipliers are matched to within a few hundred picoseconds, and the analog sum 
is formed. The resulting signal is fed to a custom waveform digitizer (WFD) with 400 
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MHz equivalent sampling rate, which provides several pulse height samples from each 
candidate electron. 

WFD data are added to the data stream only if a trigger is formed, i.e. when the 
energy associated with a pulse exceeds a pre-assigned threshold, usually taken to be 
900 MeV. When a trigger occurs, WFD samples from about 15 ns before the pulse to 
about 65 ns after the pulse are recorded. There is the possibility of two or more electron 
pulses being over-threshold in the same 80-ns time window. In that case, the length 
of the readout period is extended to include both pulses. At the earliest decay times, 
the detector signals have a large pedestal due to the 'flash,' and some of the upstream 
detectors are continuously over-threshold at early times. 

The energy and time of an electron is obtained by fitting a standard pulse shape 
to the WFD pulse using a conventional minimization. The standard pulse shape 
is established for each calorimeter. It is based on an average of the shapes of a large 
number of late-time pulses where the problems associated with overlapping pulses and 
backgrounds are greatly reduced. There are three fitting parameters, time and height 
of the pulse, and the constant pedestal, with the fits typically spanning 15 samples 
centered on the pulse. The typical time resolution of an individual electron approached 
60 ps. 

The period after each pulse is searched for any additional pulses from other 
electrons. These accidental pulses have the advantage that they do not need to be 
over the hardware threshold (~ 900 mV), but rather over the much lower (~ 250 MeV) 
minimum pulse height that can be discriminated from background by the pulse-fitting 
algorithm. A pile-up spectrum is constructed by combining the triggering pulses and 
the following accidental pulses as described later in this document. 

Zero time for a given fill is defined by the trigger pulse to the AGS kicker magnet 
that extracts the proton bunch and sends it to the pion production target. The resolution 
of the zero time needs only to be much better than the {g — 2) precession period in order 
to minimize loss of the asymmetry amplitude. 

A pulsed UV laser signal is fanned out simultaneously to all elements of the 
calorimeter stations to monitor the gain and time stabilities. For each calorimeter 
quadrant, an optical fiber carries the laser signal to a small bunch of the detector's 
scintillating fibers at the outer-radius edge. The laser pulses produce an excitation in 
the scintillators which is a good approximation to that produced by passing charged 
particles. The intensity and timing of the laser pulses are monitored with a separate 
solid-state photo-diode and a PMT, which are shielded from the beam in order to avoid 
beam-related rate changes and background which might cause shifts in the registered 
time or pulse height. The times of the laser pulses are chosen to appropriately map the 
gain and time stability during the 6 to 12 injection bunches per AGS cycle, and over 
the 10 muon lifetimes per injection. Dedicated laser runs are made 6 times per day, for 
about 20 minutes each. The average timing stability is typically found to be better than 
10 ps in any 200 /xs-interval when averaged over a number of events, with many stable 
to 5 ps. This level of timing instability contributes less than a 0.05 ppm systematic 
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error on Ua- 



3. Beam Dynamics 

The behavior of the beam in the {g — 2) storage ring directly affects the measurement of 
a^. Since the detector acceptance for decay electrons depends on the radial coordinate 
of the muon at the point where it decays, coherent radial motion of the stored beam can 
produce an amplitude modulation in the observed electron time spectrum. Resonances 
in the storage ring can cause particle losses, thus distorting the observed time spectrum, 
and must be avoided when choosing the operating parameters of the ring. Care must be 
taken in setting the frequency of coherent radial beam motion, the "coherent betatron 
oscillation" (CBO) frequency, which lies close to the second harmonic of fa = uJa/{27r). 
If /cBO is too close to 2fa the difference frequency /_ = fcBO — fa complicates the 
extraction of fa from the data, and can introduce a significant systematic error. 

A pure quadrupole electric field provides a linear restoring force in the vertical 
direction, and the combination of the (defocusing) electric field and the central magnetic 
field provides a linear restoring force in the radial direction. The {g — 2) ring is a weak 
focusing ring[43, 44, 45] with the field index 

where k is the electric quadrupole gradient. For a ring with a uniform vertical dipole 
magnetic field and a uniform quadrupole field that provides vertical focusing covering 
the full azimuth, the stored particles undergo simple harmonic motion called betatron 
oscillations, in both the radial and vertical dimensions. 

The horizontal and vertical motion are given by 

s s 
X = Xe + cos(z/^— + 4) and y = Ay cos{vy — + 5y), (24) 

Kq Kq 

where s is the arc length along the trajectory, and Rq = 7112 mm is the radius of 
the central orbit in the storage ring. The horizontal and vertical tunes are given by 
= Vl ~ n and Uy = \fri. Several n - values were used in E821 for data acquisition: 
n = 0.137, 0.142 and 0.122. The horizontal and vertical betatron frequencies are given 
by 



fx = fcVT^^ 0.929fc and /, = /cv^ ^ 0.37/c, (25) 

where fc is the cyclotron frequency and the numerical values assume that n = 0.137. 
The corresponding betatron wavelengths are Xp^ = 1.08(27ri?o) and Xjs^ = 2.7{2tiRq). It 
is important that the betatron wavelengths are not simple multiples of the circumference, 
as this minimizes the ability of ring imperfections and higher multipoles to drive 
resonances that would result in particle losses from the ring. 

The field index, ra, also determines the acceptance of the ring. The maximum 
horizontal and vertical angles of the muon momentum are given by 



^Lx = and ^^^ = ^^, (26) 
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where Xmax^Z/max = 45 mm is the radius of the storage aperture. For a betatron 
amphtude or Ay less than 45 mm, the maximum angle is reduced, as can be seen 
from the above equations. 

Resonances in the storage ring will occur if Lv^ + Mvy = N, where L, M and 
are integers, which must be avoided in choosing the operating value of the field index. 
These resonances form straight lines on the tune plane shown in Figure 17, which shows 
resonance lines up to fifth order. The operating point lies on the circle u"^ + = 1. 




Figure 17. The tunc plane, showing the three operating points used during our three 
years of running. 

For a ring with discrete quadrupoles, the focusing strength chang function 
of azimuth, and the equation of motion looks like an oscillator whose spring constant 
chang GS clS db function of azimuth s. The motion is described by 

x{s) =Xe + A^f3{s) cos{ip{s) + (5), (27) 

where (3{s) is one of the three Courant-Snyder parameters. [44] The layout of the storage 
ring is shown in Figure 10. The four-fold symmetry of the quadrupoles was chosen 
because it provided quadrupole-free regions for the kicker, traceback chambers, fiber 
monitors, and trolley garage; but the most important benefit of four-fold symmetry 
over the two-fold used at CERN[38] is that ^AnaxTAiiin = 1-03. The two-fold symmetry 

used at CERN[38] gives \J~f3^^^^^J~f3^rL = 1-15. The CERN magnetic field had significant 
non- uniformities on the outer portion of the storage region, which when combined with 
the 15% beam "breathing" from the quadrupole lattice made it much more difficult to 
determine the average magnetic field weighted by the muon distribution (Equation 20). 

The detector acceptance depends on the radial position of the muon when it 
decays, so that any coherent radial beam motion will amplitude modulate the decay 



Muon g-2: Review of Theory and Experiment 



39 



distribution. The principal frequency will be the "Coherent Betatron Frequency," 
/cBO = /c - /x = (1 - VT^)fc ^ 470 kHZ, (28) 

which is the frequency at which a single fixed detector sees the beam coherently moving 
back and forth radially. This CBO frequency is close to the second harmonic of the 
{g - 2) frequency, fa = uJa/27r ~ 228 Hz. 

An alternative way of thinking about the CBO motion is to view the ring as a 
spectrometer where the inflector exit is imaged at each successive betatron wavelength, 
A/}^. In principle, an inverted image appears at half a betatron wavelength; but the radial 
image is spoiled by the ±0.3% momentum dispersion of the ring. A given detector will 
see the beam move radially with the CBO frequency, which is also the frequency at 
which the horizontal waist precesses around the ring. Since there is no dispersion in the 
vertical dimension, the vertical waist (VW) is reformed every half wavelength A/3j^/2. A 
number of frequencies in the ring are tabulated in Table 3 

Table 3. Frequencies in the {g — 2) storage ring, assuming that the quadrupole field 
is uni form in azimuth and that n = 0.137. 



Quantity 


Expression 


Frequency 


Period 


fa 


2-Kmc^l^'^ 


0.228 MHz 


4.37 /is 


fc 


V 


6.7 MHz 


149 ns 


fx 


yi - nfc 


6.23 MHz 


160 ns 


fy 


y/nfc 


2.48 MHz 


402 ns 


/cBO 


fc fx 


0.477 MHz 


2.10 /is 


/vw 


fc 2/j, 


1.74 MHz 


0.574 /is 



The CBO frequency and its sidebands are clearly visible in the Fourier transform to 
the residuals from a fit to the five-parameter fitting function Equation 8, and are shown 
in Figure 18. The vertical waist frequency is barely visible. In 2000, the quadrupole 
voltage was set such that the CBO frequency was uncomfortably close to the second 
harmonic of fa, thus placing the difference frequency /_ = fcBO ~ fa next to /„. This 
nearby sideband forced us to work very hard to understand the CBO and how its related 
phenomena affect the value of uja obtained from fits to the data. In 2001, we carefully 
set fcBO at two different values, one well above, the other well below 2/^, which greatly 
reduced this problem. 

3.1. Monitoring the Beam Profile 

Three tools are available to us to monitor the muon distribution. Study of the beam 
de-bunching after injection yields information on the distribution of equilibrium radii in 
the storage ring. The FSDs provide information on the vertical centroid of the beam. 
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(a) (b) 



Figure 18. The Fourier transform to the residuals from a fit to the five-parameter 
fimction, showing clearly the coherent beam frequencies, (a) is from 2000, when the 
CBO frequency was close to 2a;a, and (b) shows the Fourier transform for the two 
n-values used in the 2001 run period. 

The wire chamber system and the fiber beam monitors, described above, also provide 
valuable information on the properties of the stored beam. 

The beam bunch that enters the storage ring has a time spread with a ^ 23 ns, while 
the cyclotron period is 149 ns. The momentum distribution of stored muons produces a 
corresponding distribution in radii of curvature. The distributions depend on the phase- 
space acceptance of the ring, the phase space of the beam at the injection point, and the 
kick given to the beam at injection. The narrow horizontal dimension of the beam at 
the injection point, about 18 mm, restricts the stored momentum distribution to about 
±0.3%. As the muons circle the ring, the muons at smaller radius (lower momentum) 
eventually pass those at larger radius repeatedly after multiple transits around the ring, 
and the bunch structure largely disappears after 60 /xs . This de-bunching can be seen 
in Figure 19 where the signal from a single detector is shown at two different times 
following injection. The bunched beam is seen very clearly in the left figure, with the 
149 ns cyclotron period being obvious. The slow amplitude modulation comes from the 
{g — 2) precession. By 36 /xs the beam has largely de-bunched. 

Only muons with orbits centered at the central radius have the "magic" momentum, 
so knowledge of the momentum distribution, or equivalently the distribution of 
equilibrium radii, is important in determining the correction to uja caused by the radial 
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Figure 19. The time spectrum of a single calorimeter soon after injection. The spikes 
are separated by the cyclotron period of 149 ns. 




7130 7140 7150 
Radius [mm] 



Figure 20. The distribution of equilibrium radii obtained from the beam de-bunching. 
The solid circles are from a de-bunching model fit to the data, and the dotted curve is 
obtained from a modified Fourier analysis. 



electric field used for vertical focusing. Two methods of obtaining the distribution of 
equilibrium radii from the beam debunching are employed in E821. One method uses a 
model of the time evolution of the bunch structure. A second, alternative procedure uses 
modified Fourier techniques [52]. The results from these analyses are shown in Figure 20. 
The discrete points were obtained using the model, and the dotted curve was obtained 
with the modified Fourier analysis. The two analyses agree. The measured distribution 
is used both in determining the average magnetic field seen by the muons and the radial 
electric field correction discussed below. 

The scintillating-fiber monitors show clearly the vertical and horizontal tunes as 
expected. In Figure 21, the horizontal beam centroid motion is shown, with the 
quadrupoles powered asymmetrically during scraping, and then symmetrically after 
scraping. A Fourier transform of the latter signal shows the expected frequencies, 
including the cyclotron frequency of protons stored in the ring. The traceback system 
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(b) 



Figure 21. (a) The horizontal beam centroid motion with beam scraping and without, 
using data from the scintiUating fiber hodoscopes; note the tune change between the 
two. (b) A Fourier transform of the pulse from a single horizontal fiber, which shows 
clearly the vertical waist motion, as well as the vertical tune. The presence of stored 
protons is clearly seen in this frequency spectrum. 



3.2. Corrections to uj^. Pitch and Radial Electric Field 



If the velocity is not transverse to the magnetic field, or if a muon is not at 7magic, 
the difference frequency is modified as indicated in Equation 19. Thus the measured 
frequency Ua must be corrected for the effect of a radial electric field (because of the /3 x £" 

— # — * 

term), and for the vertical pitching motion of the muons (which enters through the j3 ■ B 
term). These are the only corrections made to the uja data. We sketch the derivation, 
for E821 below[53]. For a general derivation the reader is referred to References [54, 55]. 

First we calculate the effect of the electric field, for the moment neglecting the (3 ■ B 
term. If the muon momentum is different from the magic momentum, the precession 
frequency is given by 

,Er ( 1 



(29) 



By \ a^P^^\ 
Using p = P'jm = {pm + Ap), after some algebra one finds 

i:i^ = ^ = -2^f^V (30) 

(^a By \PmJ 

Thus the effect of the radial electric field reduces the observed frequency from the simple 
frequency uJa given in Equation 14. Now 

— = (l-^)^ = (l-«)^, (31) 

Vm J^Q JM) 
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where x^ is the muon's equihbrium radius of curvature relative to the central orbit. The 
electric quadrupole field is 

nl3By 



We obtain 



E 



Auj 



KX 



Rq 



-X. 



= -2n{l - n)/3' 



UJ 



R^By 



(32) 



(33) 



so clearly the effect of muons not at the magic momentum is to lower the observed 
frequency. For a quadrupole focusing field plus a uniform magnetic field, the time 
average of x is just Xe, so the electric field correction is given by 



CE = — = -2n{l-n)f3' 

UJ 



R^By^ 



(34) 



where {xD is determined from the fast-rotation analysis (see Figure 19). The uncertainty 
on (Xg) is added in quadrature with the uncertainty in the placement of the quadrupoles 
of 5R = ±0.5 mm (±0.01 ppm), and with the uncertainty in the mean vertical position of 
the beam, ±1 mm (±0.02 ppm). For the low-n 2001 sub-period, Ce = 0.47±0.054 ppm. 



11 



Figure 22. The coordinate system of the pitching muon. The angle ip varies 
harmonically. The vertical direction is y and z is the azimuthal (beam) direction. 



The vertical betatron oscillations of the stored muons lead to l3 ■ B ^ 0. Since the 
(3 ■ B term in Equation 18 is quadratic in the components of (3, its contribution to uja will 
not generally average to zero. Thus the spin precession frequency has a small dependence 
on the betatron motion of the beam. It turns out that the only significant correction 
comes from the vertical betatron oscillation; therefore it is called the pitch correction 
(see Equation 19). As the muons undergo vertical betatron oscillations, the "pitch" 
angle between the momentum and the horizontal (see Figure 22) varies harmonically as 
ip = TpQ COS ujyt, where Uy is the vertical betatron frequency Uy = 2-71 fy, given in Equation 
25. In the approximation that all muons are at the magic 7, we set — 1/(7^ — 1) = 
in Equation 19 and obtain 



UJ„ 



m 



a^B — 



7 



>7 



{[3 ■ B)(3 



(35) 



where the prime indicates the modified frequency as it did in the discussion of the 
radial electric field given above, and Ua = —{q/m)a^B. We adopt the (rotating) 
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coordinate system shown in Figure 22, where (3 hes in the 2|/-plane, z being the 
direction of propagation, and y being vertical in the storage ring. Assuming B = yBy, 
P = z(3z + yl3y = zP cos ip + y(3 sin ip, we find 



q_ 

m ' 



a^,yBy 



7 



The small- angle approximation cosip 



ay 



and 



7 



,7 + 1, 
~ 1 and sin ip 



. 7 , 



l3yBy{zl3, + yl3y)]. (36) 
if) gives the component equations 

(37) 

(38) 



Rather than use the components given above, we can resolve u:'^ into components 
along the coordinate system defined by (3 (see Figure 22) using the standard rotation 
formula. The transverse component of uj' is given by 

(39) 



Using the small-angle expansion for cos ip 



we find 



u;_L 



2 



(40) 



Cn 




^llBy 



IpQCOS'^Uyt^ 



As can be seen from Table 3, the pitching frequency Uy is an order of magnitude larger 
than the frequency Ua, so that in one {g — 2) period uj\i oscillates more than ten times, 
thus averaging out its effect on u'^ so u'^ uj±. Thus 

2 J m 

Taking the time average yields a pitch correction 

4 nr 



(41) 



(42) 



The quantity {y^) was both 



2 4 

where we have used Equation 26 {i/jq) = n{y'^) / R^. 
determined experimentally and from simulations. For the 2001 period, Cp = 0.27 ± 
0.036 ppm, the amount the precession frequency is lowered from that given in 
Equation 20 because P ■ B 0. 

We see that both the radial electric field and the vertical pitching motion lower the 
observed frequency from the simple difference frequency Ua = {e/m)a^B, which enters 
into our determination of using Equation 16. Therefore our observed frequency must 
be increased by these corrections to obtain the measured value of the anomaly. Note 
that if ujy ^ uOa the situation is more complicated, with a resonance behavior that is 
discussed in References [54, 55]. 
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4. Analysis of the Data for uOp and uja 

In the data analysis for E821, great care was taken to insure that the results were not 
biased by previous measurements or the theoretical value expected from the standard 
model. This was achieved by a blind analysis which guaranteed that no single member 
of the collaboration could calculate the value of before the analysis was complete. 
Two frequencies, uOp, the Larmor frequency of a free proton which is proportional to 
the B field, and uJa-, the frequency that muon spin precesses relative to its momentum 
are measured. The analysis was divided into two separate efforts, oja, ujp, with no 
collaboration member permitted to work on the determination of both frequencies. 

In the first stage of each year's analysis, each independent Ua (or Up) analyzer 
presented intermediate results with his own concealed offset on Ua (or Up). Once the 
independent analyses of Ua appeared to be mutually consistent, an offset common to all 
independent Ua analyses was adopted, and a similar step was taken by the independent 
analyses of Up. The Ua offsets were kept strictly concealed, especially from the Up 
analyzers. Similarly, the Up offsets were kept strictly concealed, especially from the Ua 
analyzers. The nominal values of Ua and Up were known at best to many ppm error, 
much larger than the eventual result, and could not be guessed with any precision. No 
one person was allowed to know both offsets, and it was therefore impossible to calculate 
the value of until the offsets were publicly revealed, after all analyses were declared 
to be complete. 

For each of the four yearly data sets, 1998-2001, there were between four and five 
largely independent analyses of Ua, and two independent analyses of Up. Typically, on 
Ua there were one or two physicists conducting independent analyses in two successive 
years, and one on Up, providing continuity between the analysis of the separate data 
sets. Each of the fit parameters, and each of the potential sources of systematic error 
were studied in great detail. For the high statistics data sets, 1999, 2000 and 2001 it was 
necessary in the Ua analysis to modify the five-parameter function given in Equation 8 
to account for a number of small effects. Often different approaches were developed to 
account for a given effect, although there were common features between some of the 
analyses. 

All intermediate results for Ua were presented in terms of 3?, which is defined by 

uJa = 2TT- 0.2291MHz ■ [1 ± (3? ± A^) x 10^^] (43) 

where ±A3? is the concealed offset. Similarly, the Up analyzers maintained a constant 
offset which was strictly concealed from the rest of the collaboration. 

To obtain the value of 7?. = uja/ujp to use in Equation 16, the pitch and radial 
electric-field corrections discussed in §3.2 were added to the measured frequency Ua 
obtained from the least-squares fit to the time spectrum. Once these two corrections 
were made, the value of was obtained from Equation 16 and was published with no 
other changes. 
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The Average Magnetic Field: The ujp Analysis 

The magnetic field data consist of tliree separate sets of measurements: the cahbration 
data taken before and after each running period, maps of the magnetic field obtained 
with the NMR trolley at intervals of a few days, and the field measured by each 
of the fixed NMR probes placed outside the vacuum chamber. It was these latter 
measurements, taken concurrent with the muon spin-precession data and then tied to 
the field mapped by the trolley, which were used to determine the average magnetic 
field in the storage ring, and subsequently the value of Up to be used in Equation 16. 

4-1 .1. Calibration of the Trolley Probes The errors arising from the cross-calibration 
of the trolley probes with the plunging probes are caused both by the uncertainty in 
the relative positioning of the trolley probe and the plunging probe, and by the local 
field inhomogeneity. At this point in azimuth, trolley probes are fixed with respect 
to the frame that holds them, and to the rail system on which the trolley rides. The 
vertical and radial positions of the trolley probes with respect to the plunging probe are 
determined by applying a sextupole field and comparing the change of field measured 
by the two probes. The field shimming at the calibration location minimizes the error 
caused by the relative-position uncertainty, which in the vertical and radial directions 
has an inhomogeneity less than 0.2 ppm/cm, as shown in Fig. 23(b). The full multipole 
components at the calibration position are given in Table 4, along with the multipole 
content of the full magnetic field averaged over azimuth. For the estimated 1 mm- 
position uncertainty, the uncertainty on the relative calibration is less than 0.02 ppm. 

The absolute calibration utilizes a probe with a spherical water sample (see Figure 
14(a)) [76]. The Larmor frequency of a proton in a spherical water sample is related to 
that of the free proton through [64, 66] 

/L(sph - H20,T) = [1 - a(H20,T)] ^(free), (44) 

where (j(H20, T) is from the diamagnetic shielding of the proton in the water molecule, 
determined from [65] 

.(H,0. 34.rC, ^ 1 - ^ 

9AH) gp{H) ^p(free) 

= 25.790(14) X 10"^ (46) 

The (^-factor ratio of the proton in a spherical water sample to the electron in the 
hydrogen ground state {gj{H)) is measured to 10 parts per billion (ppb)[65]. The ratio 
of electron to proton gf-factors in hydrogen is known to 9 ppb[67]. The bound-state 
correction relating the (^-factor of the proton bound in hydrogen to the free proton are 
calculated in References [68, 69]. The temperature dependence of a is corrected for using 
da{B.20,T)/dT = 10.36(30) x 10^^/°C[70]. The free proton frequency is determined to 
an accuracy of 0.05 ppm. 

The fundamental constant A+ = yU^+//ip (see Equation 16) can be computed from 
the hyperfine structure of muonium (the iJ,~^e~ atom) [66] , or from the Zeeman splitting 
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Figure 23. Homogeneity of tlic field (a) at the calibration position and (b) for the 
azimuthal average for one trolley run during the 2000 period. In both figures, the 
contours correspond to 0.5 ppm field differences between adjacent lines. 



in muonium[39]. The latter experiment used the same cahbration probe as was used in 
our ((7 — 2) experiment; however the magnetic environments of the two experiments were 
different, so that perturbations of the probe materials on the surrounding magnetic field 
differed by a few ppb between the two experiments. 

The errors in the calibration procedure result both from the uncertainties on the 
positions of the water samples inside the trolley and the calibration probe, and from 
magnetic field inhomogeneities. The precise location of the trolley in azimuth, and 
the location of the probes within the trolley, are not known better than a few mm. 
The uncertainties in the relative calibration resulting from position uncertainties are 
0.03 ppm. Temperature and power-supply voltage dependences contribute 0.05 ppm, 
and the paramagnetism of the O2 in the air-filled trolley causes a 0.037-ppm-shift in the 
field. 



Table 4. Multipolcs at the outer edge of the storage volume (radius = 4.5 cm). The 
left-hand set arc for the plunging station where the plunging probe and the calibration 
are inserted. The right-hand set arc the multipolcs obtained by averaging over azimuth 
for a representative trolley run during the 2000 period. 
Multipole Calibration Azimuthal Averaged 

[ppm] Normal Skew Normal Skew 



Quadrupole 
Sextupole 
Octupole 
Decupole 



-0.71 
-1.24 
-0.03 
0.27 



-1.04 
-0.29 
1.06 
0.40 



0.24 
-0.53 
-0.10 
0.82 



0.29 
-1.06 
-0.15 
0.54 
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4.1.2. Mapping the Magnetic Field During a trolley run, the value of B is measured by 
each probe at approximately 6000 locations in azimuth around the ring. The magnitude 
of the field measured by the central probe is shown as a function of azimuth in Figure 
24 for one of the trolley runs. The insert shows that the fluctuations in this map that 
appear quite sharp are in fact quite smooth, and are not noise. The field maps from 
the trolley are used to construct the field profile averaged over azimuth. This contour 
plot for one of the field maps is shown in Figure 23(b). Since the storage ring has weak 
focusing, the average over azimuth is the important quantity in the analysis. Because 
NMR is only sensitive to the magnitude of B and not to its direction, the multipole 
distributions are must be determined from azimuthal magnetic field averages, where the 
field can be written as 

n=oo 

B{r, 61) = ^ r" (c„ cos nO + s„, sin n9) , (47) 

n=0 

where in practice the series is limited to 5 terms. 




Figure 24. The magnetic field measured at tlic center of the storage region vs. 
azimuthal position. Note that while the sharp fluctuations appear to be noise, when 
the scale is expanded the variations are quite smooth and represent true variations in 
the field. 



4.1.3. Tracking the Magnetic Field in Time During data-collection periods the field 
is monitored with the fixed probes. To determine how well the fixed probes permitted 
us to monitor the field felt by the muons, the measured field, and that predicted by 
the fixed probes is compared for each trolley run. The results of this analysis for the 
2001 running period is shown in Figure 25. The rms distribution of these differences is 
0.10 ppm. 

4.1.4. Determination of the Average Magnetic Field: uOp The value of Up entering 
into the determination of is the field profile weighted by the muon distribution. The 
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Figure 25. The difference between the average magnetic field measured by the trolley 
and that inferred from tracking the magnetic field with the fixed probes between trolley 
maps. The vertical lines show when the magnet was powered down and then back up. 
After each powering of the magnet, the field does not exactly come back to its previous 
value, so that only trolley runs taken between magnet powerings can be compared 
directly. 



multipoles of the field, Equation 47, are folded witli tlie muon distribution 

M(r, 6) = ^[7m(^) cosmO + ^^(r) sin m6'], (48) 

to produce the average field, 

(5)^_dist = / M(r, 9)B{r, 9)rdrd9, (49) 

where the moments in the muon distribution couple moment-by-moment to the 
multipoles of B. Computing {B) is greatly simplified if the field is quite uniform (with 
small higher multipoles), and the muons are stored in a circular aperture, thus reducing 
the higher moments of M(r, 9). This worked quite well in E821, and the uncertainty on 
{B) weighted by the muon distribution was ±0.03 ppm. 

The weighted average was determined both by a tracking calculation that used a 
field map and calculated the field seen by each muon, and also by using the quadrupole 
component of the field and the beam center determined from a fast-rotation analysis to 
determine the average field. These two agreed extremely well, vindicating the choice of 
a circular aperture and the ±1 ppm specification on the field uniformity, that were set 
in the design stage of the experiment. [26] 

4-1.5. Summary of the Magnetic Field Analysis The limitations on our knowledge 
of the magnetic field come from measurement issues, not statistics, so in E821 the 
systematic errors from each of these sources had to be evaluated and understood. The 
results and errors are summarized in Table 5. 
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Table 5. Systematic errors for the magnetic field for the different run periods. 
^Higher multipoles, trolley temperature and its power-supply voltage response, and 
eddy currents from the kicker. 



Source of errors 


1999 


2000 


2001 




fppml 


fppml 


fppml 


Absolute calibration of standard probe 


0.05 


0.05 


0.05 


Calibration of trolley probes 


0.20 


0.15 


0.09 


Trolley measurements of Bq 


0.10 


0.10 


0.05 


Interpolation with fixed probes 


0.15 


0.10 


0.07 


Uncertainty from muon distribution 


0.12 


0.03 


0.03 


Inflector fringe field uncertainty 


0.20 






Otherf 


0.15 


0.10 


0.10 


Total systematic error on ujp 


0.4 


0.24 


0.17 


Muon-averaged field [Hz]: uj^jli^ 


61791256 


61791595 


61791400 



The Muon Spin Precession Frequency: The Ua Analysis 
To obtain the muon spin precession frequency uJa given in Equation 14, 
^ qB 

= -af, — , (50) 
m 

which is observed as an oscillation of the number of detected electrons with time 

N{t, Eth) = N^{Eth)e-"^^[l + A{Eth) cos{uJat + 0(^th))], (51) 
it is necessary to: 

• Modify the five-parameter function above to include small effects such as the 
coherent betatron oscillations (CBO), pulse pile- up, muon losses, and gain changes, 
without adding so many free parameters that the statistical power for determining 
oja is compromised. 

• Obtain an acceptable per degree of freedom in all fits, i.e. consistent with 1, 
where ct(x|) = ^2/NDF. 

• Insure that the fit parameters are stable independent of the starting time of the 
least-square fit. This was found to be a very reliable means of testing the stability 
of fit parameters as a function of the time after injection. 

In general, fits are made to the data out to about 640 /is, about 10 muon lifetimes. 

4-2.1. Distribution of Decay Electrons Decay electrons with the highest laboratory 
energies, typically E > 1.8 GeV, are used in the analysis for Ua, as discussed in § 1.1. 
In the (excellent) approximations that (3 ■ B = and the effect of the electric field on 
the spin is small, the average spin direction of the muon ensemble (i.e., the polarization 
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vector) precesses, relative to the momentum vector, in the plane perpendicular to the 
magnetic field, B = Byf), according to 

S = (S_L sin {Uat + (f>)x + Syij + S±_ cos {ijJat + (j))z). (52) 

The unit vectors £, and z are directed along the radial, vertical and azimuthal 
directions respectively. The (constant) components of the spin parallel and 
perpendicular to the B-field are Sy «1 and with ^ {s\ + Sy) = 1. The polarization 
vector precesses in the plane perpendicular to the magnetic field at a rate independent of 
the ratio Sy/s±. Since = {g — 2)/2 > 0, the spin vector rotates in the same plane, but 
slightly faster than the momentum vector. Note that the present experiment is, apart 
from small detector acceptance effects, insensitive to whether the spin vector rotates 
faster or slower than the momentum vector rotation, and therefore it is insensitive to 
the sign of a^. There are small geometric acceptance effects in the detectors which 
demonstrate that our result is consistent with > 0. 

The value for Ua is determined from the data using a least-square minimization 
fit to the time spectrum of electron decays, = J2i i^i ~ ^ / N {ti) , where the iVj 
are the data points, and N{ti) is the fitting function. The statistical uncertainty, in 
the limit where data are taken over an infinite number of muon lifetimes, is given by 
Equation 11. The statistical figure of merit is Fm = A\fN which reaches a maximum 
at about y = 0.8, ot E ^ 2.6 GeV/c (see Figure 2). If all electrons are taken above 
some minimum energy threshold, FM{Eth) reaches a maximum at about y = 0.6, or 
Ethresh = 1.8 GeV/c (Figure 3). 

The spectra to be fit are in the form of histograms of the number of electrons 
detected versus time (see Figure 26), which in the ideal case follow the five-parameter 
distribution function. Equation 8. While this is a fairly good approximation for the 
E821 data sets, small modifications, due mainly to detector acceptance effects, must be 
made to the five-parameter function to obtain acceptable fits to the data. The most 
important of these effects are described in the next section. 

The five-parameter function has an important, well-known invariance property. A 
sum of arbitrary time spectra, each obeying the five-parameter distribution and having 
the same A and u, but different values for A'o, A, and 0, also has the five-parameter 
functional form with the same values for A and u. That is, 

^5,e~^(*'^*«)(l+AiCos {uj{ti - tio) + (pi)) = 5e-^*(l+Acos(cjt + 0)).(53) 

i 

This invariance property has significant implications for the way in which data are 
handled in the analysis. The final histogram of electrons versus time is constructed from 
a sum over the ensemble of the time spectra produced in individual spills. It extends in 
time from less than a few tens of microseconds after injection out to 640/is, a period of 
about 10 muon lifetimes. To a very good approximation, the spectrum from each spill 
follows the five-parameter probability distribution. From the invariance property, the 
t=0 points and the gains from one spill to the next do not need to be precisely aligned. 

Pulse shape and gain stabilities are monitored primarily using the electron data 
themselves rather than laser pulses, or some other external source of pulses. The electron 
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Figure 26. Histogram of the total number of electrons above 1.8 GeV versus time 
(modulo 100 jjL s) from the 2001 /i^ data set. The bin size is the cyclotron period, 
K, 149.2 ns, and the total number of electrons is 3.6 billion. 
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Figure 27. Typical calorimeter energy distribution, with an endpoint fit 
superimposed. The inset shows the full range of reconstructed energies, from 0.3 to 
3.5 GeV. 



times and energies are given by fits to standard pulse shapes, which are are estabhshed 
for each detector by taking an average over many pulses at late times. The variations in 
pulse shapes in all detectors are found to be sufficiently small as a function of energy and 
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decay time, and contribute negligibly to the uncertainty in uja- In order to monitor the 
gains, the energy distributions integrated over one spin precession period and corrected 
for pile-up are collected at various times relative to injection. The high-energy portion 
of the energy distribution is well-described by a straight line between the energy points 
at heights of 20% and 80% of the plateau in the spectrum (see Figure 27). The position 
of the X-axis intercept is taken to be the endpoint energy, 3.1 GeV. It is found that the 
energy stability of the detectors on the 'quiet' side of the ring stabilize earlier in the 
spill than those on the noisy side of the ring. Some of the gain shift is due to the PMT 
gating operation. Since the noisy detectors are gated on later in the spill than the quiet 
ones, their gains tend to stabilize later. The starting times for the detectors are chosen 
so that most of their gains are calibrated to better than 0.2%. On the quiet side of the 
ring, data fitting can begin as early as a couple of microseconds after injection; however 
it is necessary to delay at least until the beam-scraping process is completed. For the 
noisiest detectors, just downstream of the injection point, the start of fitting may be 
delayed to 30 /xs or more. The time histograms are then accumulated after applying the 
gain correction to the energy of each electron. An uncertainty in the gain stability on 
average over a fill affects A^, A, 0, and A to a small extent. The result is a systematic 
error on uja on the order of 0.1 ppm. 

The cyclical motion of the bunched beam around the ring (fast-rotation structure, 
§3.1) results in an additional multiplicative oscillation superimposed on the usual five- 
parameter spectrum, which is clearly visible in Figure 19. The amplitude of the 
modulation, which we refer to as the 'fast-rotation' structure, dies away with a lifetime 
~ 26/is as the muon bunch spreads out around the ring (see § 3.1). In order to filter the 
fast-rotation structure from the time spectrum, two steps are taken. First, the t = 
point for electron data from each spill is offset by a random time chosen uniformly 
over one cyclotron period of 149 ns. Second, the data are formed into histograms with 
a bin width equal to the cyclotron period. The fast-rotation structure is also partially 
washed out when data are summed over all detectors, since the fast-rotation phase varies 
uniformly from to 27r in going from one detector to the next around the ring. These 
measures suppress the fast-rotation structure from the spectrum by better than a factor 
of 100. It is confirmed in extensive simulations that these fast rotation filters have no 
significant effect on the derived value for uja- We note that the accidental overlap of 
pulses ('pile-up') is enhanced by the fast-rotation structure and is accounted for in the 
process of pile-up reconstruction, as described later in this section. 

The five-parameter functional shape is unaffected by the size of the bin width in 
the histogram. The number of counts in each bin of the time histogram is equal to 
the integral of the five-parameter function over the bin width. The integrated five- 
parameter function has the same A and uj as the differential five-parameter function by 
the invariance property. However, we note that unduly large values for the bin width 
will reduce the asymmetry and perhaps introduce undesirable correlations between uja 
and the other parameters of the fit. A bin width equal to the cyclotron period of about 
149 ns is chosen in most of the analyses to help suppress the fast-rotation structure. 
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This is narrow enough for binning effects to be minimaL 

4-2.2. Modification of the Five-Parameter Fitting Function In order to successfully 
describe the functional form of the spectrum of detected electrons versus time, the five- 
parameter function must be modified to include additional small effects of detector 
acceptance, muon losses, electron drift time from the point of muon decay to the 
point of detection, pile-up, etc. Given the desire to maintain as nearly as possible 
the invariance property of the five-parameter function, and the fact that some of the 
necessary corrections to the function turn out to be imprecisely known, every effort 
has been made to design the apparatus so as to minimize the needed modifications. 
Fortunately, most of the necessary modifications to the five-parameter function lead 
to additional parameters having little correlation with Ua and therefore contribute 
minimally to the systematic error. The following discussion will concentrate on those 
modifications that contribute most to the systematic error on Ua- 

Two separate decay electrons arriving at one calorimeter with a time separation 
less than the time resolution of the calorimeters (~ 3 to 5 ns, depending on the pulse 
heights) can be mistaken for a single pulse (pile-up). The registered energy and time in 
this case will be incorrect, and if not properly accounted for, it can be a source of a shift 
in uja. The registered energy and time of a pileup pulse is approximately the sum of 
the energies and the energy-weighted average time respectively of the two pulses. The 
phase of the {g — 2) precession depends on energy, and since the pile- up pulse has an 
incorrectly registered energy, it has an incorrect phase. Such a phase shift would not be 
a problem if pile-up were the same at early and late decay times. However, the rate of 
pile-up is approximately proportional to the square of the total decay rate, and therefore 
the ratio of its rate to the normal data rate varies as ~ e~*^"'^. The ear ly-to- late change 
in relative rates leads to an early-to-late shift in the pile- up phase (ppuit)- When the 
oscillating term cos{uJat + (ppuit) + <f)) is fit to cos{uJat + 0), a shift in the value for Ua 
can result. 

To correct for pile-up-induced phase shifts, an average pile-up spectrum is 
constructed and then subtracted from the data spectrum. It is constructed as follows. 
When a pulse exceeds a pre-assigned energy threshold, it triggers the readout of a regular 
electron pulse, which consists of about 80 ns of WFD information with a pulse height 
reading every 2.5 ns in that time interval. This time segment is large enough to include 
the trigger pulse (S'r(i?,t)) plus a significant period beyond it, which usually contains 
only pedestal. Occasionally, a second ('shadow', Ss{E,t)) pulse from another electron 
falls in the pedestal region. If it falls into a pre-selected time interval after the trigger 
pulse, it is used to construct the pileup spectrum by combining the 5*^ and Ss pulses 
into a single pile-up pulse, D. The width of the pre-selected time interval is chosen to be 
equal to the minimum time separation between pulses which the detectors are capable 
of resolving (typically around 3 ns, however this depends slightly on the relative energies 
in the two peaks). The fixed delay from the trigger pulse and the window is subtracted 
from the time of Ss and then the pulse D is formed from the sum of the WFD samples 



Muon g-2: Review of Theory and Experiment 



55 



of St and the Ss- The pile-up spectrum is formed from P{E,t) = D — St — Ss- By- 
construction, the rate of shadow pulses per unit time is properly normalized to the 
trigger rate to a good approximation. The single pulses are explicitly subtracted when 
forming P{E, t) because two single pulses are lost for every pile-up pulse produced. To 
properly account for the impact on the pile-up due to the rapid variation in data rate 
produced by the fast-rotation structure, the time difference between the trigger pulse 
and the shadow pulse's time window is kept small compared to the fast-rotation period 
of 149 ns. 

The resulting pile-up spectrum is then subtracted from the total spectrum, on 
average eliminating the pile-up. As Figure 28 demonstrates, the constructed pile-up 
spectrum is in excellent agreement with the data above the endpoint energy of 3.1 
GeV, where only pile-up events can occur. Below about 2.5 GeV, the pile-up spectrum 
becomes negative because the number of single pulses, St and 5*5 lost exceeds the 
number of pile-up D pulses gained. 

The magnitude of the pile-up is reconstructed with an estimated uncertainty of 
about ±8%. There is also an uncertainty in the phase of the pile-up. Further, the 
pulse-fitting procedure only recognizes shadow pulses above about 250 MeV. Based on 
simulations, the error due to 'unseen pile-up' (pulses below 250 MeV) was about 0.03 
ppm in the 2001 data set. The systematic uncertainties were 0.014 ppm and 0.028 
ppm from uncertainties in the pile-up amplitude and phase respectively, and a total 
contribution of 0.08 ppm came from pile-up effects in the 2001 data set. 




Figure 28. Absolute value of the constructed pile-up spectrum (solid line) and actual 
data spectrum (dashed line). The two curves agree well at energies well above the 
electron endpoint at 3.1 GeV, as expected. The constructed spectrum is negative below 
2.5 GeV because the number of single electrons lost to pile- up exceeds the number of 
pile-up pulses. 

The coherent betatron oscillations, or oscillations in the average position and 
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width of the stored beam, (§3) cause unwanted oscillations in the muon decay time 
spectrum. These effects are generically referred to as CBO. Some of the important 
CBO frequencies are given in Table 3. They necessitate small modifications to the five- 
parameter functional form of the spectrum, and, like pile-up, can cause a shift in the 
derived value of Ua if they are not properly accounted for in the analysis. 

In particular, the functional form of the distribution of detector acceptance versus 
energy changes as the beam oscillates, affecting the number and average energy of the 
detected electrons. The result is that additional parameters need to be added to the 
five-parameter function to account for CBO oscillations. The effect of vertical betatron 
oscillations is small and dies away much faster than the horizontal oscillations, so that 
it can be neglected at the fit start times used in most of the analyses. 

The horizontal CBO affects the spectra mainly in two ways: it causes oscillations in 
the number of particles because of an oscillation in the acceptance of the detectors, and 
it causes oscillations in the average detected energy. For a time spectrum constructed 
from the decay electrons in a given energy band, oscillation in the parameter N is due 
primarily to the oscillation in the the detector acceptance. Oscillations induced in A 
and 0, on the other hand, depend primarily on the oscillation in the average energy. In 
either case, each of the parameters A^, A and acquire small CBO-induced oscillations 
of the general form 

Pi = Pio[l + B,e-^^^° cos {ucBot + 9,)]. (54) 

The presence of the CBO introduces fcBOi its harmonics, and the sum and 
difference frequencies associated with beating between fcBO and fa = oJa/^Tx ^ 229.1 
kHz. If the CBO effects are not included in the fitting function, it will pull the value of 
uja in the fit by an amount related to how close fcBO is to the second harmonic of /„, 
introducing a serious systematic error (see Table 3). This effect is shown qualitatively 
in Figure 29. 




Figure 29. The relative shift in the value obtained for Ua as a function of the CBO 
frequency, when the CBO effects are neglected in the fitting function. The vertical line 
is at 2/a, and the operating points for the different data collection periods is indicated 
on the curve. 
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The problems posed by the CBO in the fitting procedure were solved in a variety 
of ways in the many independent analyses. All analyses took advantage of the fact 
that the CBO phase varied fairly uniformly from to 2-7r around the ring in going from 
one detector to the next; the CBO oscillations should tend to cancel when data from all 
detectors are summed together. In the 2001 data set, where the field index n was chosen 
to move the difference frequency /_ = /cbo ~ fa well away from fa, one of the analyses 
relied only on the factor-of-nine reduction in the CBO amplitude in the sum of data 
over all detectors, and no CBO-induced modification to the five-parameter function was 
necessary. 

The cancellation of CBO from all detectors combined would be perfect if all the 
detector acceptances were identical, or even if opposite pairs of detectors at 180° in 
the ring were identical. Imperfect cancellation is due to the reduced performance of 
some of the detectors and to slight asymmetries in the storage-ring geometry. This was 
especially true of detector number 20, where there were modifications to the vacuum 
chamber and whose position was displaced to accommodate the traceback chambers. 
The 180° symmetry is broken for detectors near the kicker because the electrons pass 
through the kicker plates. Also the fit start times for detectors near the injection point 
are inevitably later than for detectors on the other side of the ring because of the 
presence of the injection flash. 

In addition to relying on the partial CBO cancellation around the ring, all of 
the other analysis approaches use a modified function in which all parameters except 
uja and A oscillate according to equation 54. In fits to the time spectra, the CBO 
parameters ojcbo and \cbo ~ 100/xs (the frequency and lifetime of the CBO oscillations, 
respectively), are typically held fixed to values determined in separate studies. They 
were established in fits to time spectra formed with independent data from the FSDs 
and calorimeters, in which the amplitude of CBO modulation is enhanced by aligning 
the CBO oscillation phases of the individual detectors and then adding all the spectra 
together. 

Using these fixed values for ojcbo and Acso, fits to the regular spectra are made 
with PjO; Bi and (f)i in Equation 54 (one set of these three for each of A^, A, </>) as free 
parameters. With the addition of Ua and A as parameters, one obtains a total of 11 free 
parameters in the global fit to data. In practice, in some of the analyses, the smaller 
CBO parameters are held fixed or eliminated altogether. 

Another correction to the fitting function is required to account for muons being 
lost from the storage volume before they have had a chance to decay. Such losses lead to 
a distortion of the spectrum which can result in an incorrect fitted value for the lifetime 
and a poor value for the reduced of the fit. The correlation between the lifetime 
and the precession frequency is quite small, so that the loss in number of counts does 
not really affect the value of the precession frequency. The major concern is that the 
average spin phase of muons lost might be different from that of the muons which remain 
stored. If this were the case, the average phase of the stored muons would shift as a 
function of time, leading directly to an error in the fitted value of uJa- This uncertainty 
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in the phase forces an inflation in the systematic error due to muon losses. Muon losses 
from the ring are believed to be induced by the coupling between the higher moments 
of the muon and the field distributions, driving resonances that result in an occasional 
muon striking a collimator and leaving the ring. The losses are as large as 1% per muon 
lifetime at early decay times, but typically settle down to less than 0.1% 200 /is after 
injection. Such losses require the modification Nq Nq{1 — Aioss'^ioss(^)) in the fitting 
function, with Ai^ss being an additional parameter in the fit. The function riioss, which 
represents the time distribution of lost muons, is obtained from the muon loss monitors. 

An alternative analysis method to determine uja utilizes the so-called 'ratio method.' 
Each electron event is randomly placed with equal probability into one of four time 
histograms, Ni — N4, each looking like the usual time spectrum. Figure 26. A spectrum 
based on the ratio of combinations of the histograms is formed: 

For a pure five-parameter distribution, keeping only the important large terms, this 
reduces to 

r{t) = A cos {ooat + 0) + (56) 

where Ta = ^Ti/ua is an estimate (~ 10 ppm is easily good enough) of the spin 
precession period, and the small constant offset produced by the exponential decay 
is yq{-^Y = 0.000287. Construction of independent histograms and A'4 simplifies 
the estimates of the statistical uncertainties in the fitted parameters. 

The advantage of fitting r(t), as opposed to the standard technique of fitting N{t), 
is that any slowly varying (e.g. with a period much larger than r^) multiplicative 
modulation of the five-parameter function, which includes the exponential decay of the 
muon itself, largely cancels in this ratio. Other examples of slow modulation are muon 
losses and small shifts in PMT gains due to the high rates encountered at early decay 
times. There are only 3 parameters. Equation 56, compared to 5, Equation 8, in the 
regular spectrum. Note that faster-varying effects, such as the CBO, will not cancel in 
the ratio and must be handled in ways similar to the standard analyses. One of the ratio 
analyses of the 2001 data set used a fitting function formed from the ratio of functions 
hi, consisting of the five-parameter function modified to include parameters to correct 
for acceptance effects such as the CBO: 

^ hijt + iTg) + /i2(t - |r,) - hsjt) - hjjt) 

- h,{t + ir,) + h2it - ir,) + h,it) + h,it) ^ ' 

The systematic errors for three yearly data sets, 1999 and 2000 for and 2001 for 
fi~ , are given in Table 6. 
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Table 6. Systematic errors for uJa in the 1999, 2000 and 2001 data periods. In 
2001, systematic errors for the AGS background, timing shifts, £'-field and vertical 
oscillations, beam de-bunching/randomization, binning and fitting procedure together 
equaled 0.11 ppm and this is indicated by X in the table. 



"^syst 


1999 


2000 


2001 




(ppm) 


(ppm) 


(ppm) 


Pile-up 


0.13 


0.13 


0.08 


AGS Background 


0.10 


0.01 


t 


Lost Muons 


0.10 


0.10 


0.09 


Timing Shifts 


0.10 


0.02 


X 


E-field and Pitch 


0.08 


0.03 


X 


Fitting/Binning 


0.07 


0.06 


X 


CBO 


0.05 


0.21 


0.07 


Gain Changes 


0.02 


0.13 


0.12 


Total for uOa 


0.3 


0.31 


0.21 



5. The Permanent Electric Dipole Moment of the Muon 

If the muon were to possesses a permanent electric dipole moment (see Equation 2), an 
extra term would be added to the spin precession equation (Equation 20): 



^EDM — 77 

m I 



0xB + -iE-^0-E)0) 
c 7 + 1 



(5^ 



where the dimensionless constant rj is proportional to the muon electric dipole moment, 
d = V2^^- The dominant (3 x B term is directed radially in the storage ring, transverse 
to UJa- The total precession vector, uj = uja + ^edm, is tipped from the vertical direction 
by the angle 5 = tan~^ Equation 52, with the simplification that initially Sy = 0, 
s_i_ = 1 and = 0, is modified to 

s = (cos 6 sin Uat x + sinS sin Uat y + cos Uat z) . (59) 

The tipping of the precession plane produces an oscillation in the average vertical 
component of the spin which, because of the correlation between the spin and electron 
momentum directions, in turn causes oscillation in the average vertical component of 
the electron momentum. The average vertical position of the electrons, at the entrance 
face of the calorimeters, oscillates at angular frequency 90 degrees out of phase with 
the number oscillation depicted in Figure 26, with an amplitude proportional to the 
EDM. Additionally, the magnitude of of the precession frequency is increased by the 
EDM, 



A measure (or limit) of the muon EDM, independent of the value of a^, will be 
produced from an on-going analysis of the vertical electron motion using the FSD, PSD 
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and traceback data from E821. Furthermore, a difference between tlie experimental and 
standard model values of could be generated by an electric dipole moment. 

No intrinsic electric dipole moment (EDM) has ever been experimentally detected 
in the muon or in any other elementary particle or atom. The existence of a permanent 
(EDM) in an elementary particle would violate both T (time reversal) and P (parity) 
symmetries. This is in contrast to the magnetic dipole moment (MDM), which is allowed 
by both of these symmetries. A non-vanishing EDM (or MDM) is consistent with C, 
charge conjugation symmetry, and with the combined symmetry, CPT, provided that 
the magnitudes of the EDM (or MDM) for a particle and its anti-particle be equal in 
magnitude but opposite in sign. No violation of CPT symmetry has ever been observed, 
and it is strictly invariant in the Standard Model. Assuming CPT invariance, the T 
and P violations of a non- vanishing EDM imply a violation of CP and CT, respectively. 
While P violation in the weak interactions is maximal and has been observed in many 
reactions, CP violation has only been observed in decays of neutral K and B mesons, 
with very small amplitudes. The violation of T invariance has only been observed in the 
decays of neutral K mesons [78]. An EDM large enough to be detected would require 
new physics beyond the standard model. 

It is widely believed that CP violation is one of the ingredients needed to explain the 
baryon asymmetry of the universe, however the type and size of CP violation which has 
been experimentally observed to date is much too small to account for it. Searches for 
EDMs are the subject of many past, current and future experiments since its detection 
would herald the presence of new sources of CP violation beyond the Standard Model. 
Indeed many proposed extensions to the Standard Model predict or at least allow the 
possibility of EDMs and the failure to observe them has historically served as a powerful 
constraint on these models. 

The current experimental limit on the muon EDM, a by-product of the third CERN 
{g — 2) experiment, is < 1 x 10^^^ e-cm (90% C.L.)[63]. This is based on the non- 
observation of an oscillation in the number of electrons above compared to below the 
calorimeter mid-planes at the precession frequency, fa, 90 degrees out of phase with the 
number oscillation. 

A new muon EDM limit will be set using E821 data, by searching for an oscillation 
in the average vertical positions of electrons from the FSD and PSD data. The new 
data have the advantage over the CERN data of containing information on the shape of 
the vertical distribution of electrons, not just the number above or below the detector 
mid-plane, and the sizes of the data samples are significantly larger. Analyses of these 
data are in progress. 

A major systematic error in the EDM measurements by CERN, and by E821 using 
the PSD and FSD information, arises from any misalignment between the vertical center 
of the stored muons and the vertical positions of the detectors. As the spin precesses, 
the vertical distribution of electrons entering the calorimeters changes. In the case of the 
H~ decays, high energy (in the MRF) electrons are more likely to be emitted opposite 
rather than along the muon spin direction. When the spin has an inward (outward) 
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radial component in the storage ring, the average electron will have a positive (negative) 
radial component of laboratory momentum, and will have to travel a greater (lesser) 
distance to get to the detectors. When the spin points inward, more electrons are 
emitted outward, and on average the electrons travel further and spread out more in 
the vertical direction before getting to the detectors. The change in vertical distribution 
combined with a misalignment between the beam and the detectors, inevitably leads to 
a vertical oscillation which appears as a false EDM signal. This is the major limitation 
to this approach for measuring the EDM, and substantial improvements in the EDM 
limit in the future will require a dedicated experiment to circumvent this problem[79]. 

In E821, the traceback chambers are able to measure oscillation of the vertical 
angles of the electron trajectories as well as average position, providing a limit on the 
muon EDM which is complementary to that obtained from the FSD and PSD data. The 
angle information from the traceback data is less susceptible to the systematic error from 
beam misalignment than up-down data, however the statistical sample is smaller than 
the FSD or PSD data. It is expected that the final EDM result from the E821 the FSD, 
PSD, and traceback data will improve upon the present limit on the muon EDM by a 
factor of four to five, to the level of ~ few x 10~^^ e-cm. 

In the process of deducing the experimental value of (see equation 21) from the 
measured values of uja and cUp, it is assumed that the muon EDM is negligibly small and 
can be neglected. It is found under this assumption and discussed in a later section of 
this manuscript, that the experimental value of the anomaly differs from the standard 
model prediction (equation 162), 

Aa^ = 295(88) x 10"^^ (61) 

The extreme view could be taken that this apparent difference, Aa^, is not due to 
a shift in the magnetic anomaly, but rather is due entirely to a shift of the precession 
frequency, cj^, caused by a non-vanishing EDM. From Equation 60 the EDM would have 
to have the value = 2.4(0.4) x 10~^^ e-cm. This is a a factor of ~ 10® larger than the 
current limit on the electron EDM. Such a large muon EDM is not predicted by even 
the most speculative models outside of the standard model, but cannot be excluded by 
the previous CERN experimental limit on the EDM, and also may not be definitively 
excluded by the eventual E821 experimental result from the PSD, FSD and traceback 
data. Of course, if the change in the precession frequency were due to an EDM rather 
than a shift in the anomaly, then this would be also be a very interesting indication of 
new physics [80]. 

6. Results and Summary of E821 

The values obtained in E821 for are given in Table 1. However, the experiment 
measures TZ = Ua/^p, not directly, where the tilde over Ua means that the pitch and 
radial electric field corrections have been included (see §3.2). The fundamental constant 
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A+ = fi^+/ fip (see Equation 16) connects the two quantities. The values obtained for TZ 
are given in Table 7. 

Table 7. The frequencies uJa and Wp obtained from the three major data-collection 
periods. The radial electric field and pitch corrections applied to the Wa values are 
given in the second column. Total uncertainties for each quantity are shown. The 
right-hand column gives the values of TZ, where the tilde indicates the muon spin 
precession frequency corrected for the radial electric field and the pitching motion. 
The error on the average includes correlations between the systematic uncertainties of 
the three measurement periods. 



Period 



Ua/i2n) [Hz] 



£'/pitch [ppm] ujp/{2n) [Hz] 



1999 (/i- 

2000 {[T 

2001 {^i 



229 072.8(3) 
229 074.11(16) 
229 073.59(16) 



+0.81(8) 
+0.76(3) 
+0.77(6) 



61791256(25) 
61791595(15) 
61791400(11) 



0.003 7072041(51) 
0.003 707205 0(2 5) 
0.003 707208 3(2 6) 



Average 



0.003 707206 3(2 0) 



The results are 



and 



7^^- = 0.003 707208 3(2 6) 



7^ + = 0.003 707204 8(2 5). 



(62) 



(63) 



The CPT theorem predicts that the magnitudes of 7^^+ and TZ^- should be equal. The 
difference is 



A7^ = lZ^,r- - TZ^+ = (3.5 ± 3.4) x 10" 



(64) 



Note that it is the quantity TZ that must be compared for a CPT test, rather than 
since the quantity A+ which connects them is derived from measurements of the 



hyperfine structure of the jj,~^e~ atom (see Equation 16). 

Using the latest value A+ = /i^+/yUp = 3.183 345 39(10) [39] gives the values for 
shown in Table 1. Assuming CPT invariance, we combine all measurements of and 
a~ to obtain the 'world average' [26] 

a/,(Expt) = 11659 208.0(6.3) X 10"^° (0.54ppm). (65) 

The values of obtained in E821 are shown in Figure 30. The final combined 
value of represents an improvement in precision of a factor of 13.5 over the CERN 
experiments. The final error of 0.54 ppm consists of a 0.46 ppm statistical component 
and a 0.28 systematic component. 
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Figure 30. Results for the E821 individual measurements of by running year, 
together with the final average. 



7. The Theory of the Muon g-2 

As discussed in the introduction, the gf-factor of the muon is the quantity which relates 

its spin s to its magnetic moment /I in appropriate units: 

q 

^^^^2m~^' and gjj^=2{l + a^) . (66) 

Dirac 

In the Dirac theory of a charged spin-1/2 particle, g = 2. Quantum Electrodynamics 
(QED) predicts deviations from the Dirac exact value, because in the presence of an 
external magnetic field the muon can emit and re-absorb virtual photons. The correction 
Qfj, to the Dirac prediction is called the anomalous magnetic moment. As we have seen 
in the previous section, it is a quantity directly accessible to experiment. 

In this section, we shall present a review of the various contributions to in the 
Standard Model, with special emphasis on the evaluations of the hadronic contributions. 

7.1. The QED Contributions 

In QED of photons and leptons alone, the Feynman diagrams which contribute to at 
a given order in the perturbation theory expansion (powers of ^), can be divided into 
four classes: 

7.1.1. Diagrams with Virtual Photons and Muon Loops. 

Examples are the lowest-order contribution in Figure 31 and the two- loop contributions 
in Figure 32. In full generality, this class of diagrams consists of those with virtual 
photons only (wavy black lines), and those with virtual photons and internal fermion 
loops (solid blue loops) restricted to be of the same fiavor as the external line (solid 
blue line) in an external magnetic field (X in the diagrams). Since is a dimensionless 
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quantity, and there is only one type of fermion mass in these graphs-with no other mass 
scales- these contributions are purely numerical, and they are the same for the three 
charged leptons: / = e,/i,r. Indeed, a] from Figure 31 is the celebrated Schwinger 
result [16], mentioned previously. 




Figure 31. Lowest-order QED contribution. The solid blue line represents the muon 
in an external magnetic field (X in the figure). The wavy black line represents the 
virtual photon. 

while af^ from the seven diagrams in Figure 32 (the factor of 2 in a diagram corresponds 
to the contribution from the mzrror diagram) gives the result[81, 82] 




At the three-loop level there are 72 Feynman diagrams of this type. Only a 
few representative examples are shown in Figure 33. Quite remarkably, their total 
contribution is also known analytically (see Reference [83] and references therein). They 
bring in transcendental numbers like C(3) = 1.202 0569 ■ ■ ■, the Riemann zeta-function 
of argument 3 which already appears at the two-loop level in Equation (68), as well as 
transcendent als of higher complexity: * 



.(6) 



28259 17101 



5184 
100 



810 



2 298 2, 



Li4(l/2) + - 



TT 



139 



239 



18 ^^^^ 2160 



TT 



* There is in fact an interesting relationship between the appearance of these transcendentals in 
pcrturbativc quantum field theory and mathematical structures like knot theory and non-commutative 
geometry, which is under active study (see, e.g., References [84, 85] and references therein). 
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Figure 32. Two-loop QED Fcynman diagrams with the same lepton flavor. 









Figure 33. A few Feynman diagrams of the three-loop type. In this class the flavor 
of the internal fermion loops is the same as the external fermion. 



|'^«3)-fc(5)}e) 



1.181 241 456, 



a 



(70) 



where Li4(l/2) = 0.517 479 ■ ■ -, is a particular case of the polylogarithm-function (see, 
e.g. References [86, 87]) which often appears in loop calculations in perturbation theory, 
in the form of the integral representation: 
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Lifc(x) = 1^-^ |ln^-2tln(l - xt) , k>2 (71) 

OO / 1 N fc 

= E - N<i> (72) 



n=l 



while 

OO -| 

C(fc)=Li,(l) = ^-; A;>2. (73) 

At the four-loop level, there are 891 Feynman diagrams of this type. Some of them 
are already known analytically, but in general one has to resort to numerical methods for 
a complete evaluation. This impressive calculation, which requires many technical skills 
(see the chapter Theory of the anomalous magnetic moment of the electron-Numerical 
Approach in Reference [88] for an overall review), is under constant updating due to 
advances in computing technology. The most recent published value for the whole four- 
loop contribution, from fermions with the same flavor, gives the result (see Reference [89] 
for details and references therein) 

f8)_ 



-1.728 3(35) (^-j , (74) 

where the error is due to the present numerical uncertainties. 

Notice the alternating sign of the results from the contributions of one loop to four 
loops, a simple feature which is not yet a priori understood . Also, the fact that the sizes 
of the (^^^ coefficients for n = 1,2,3,4 remain rather small is an interesting feature, 
allowing one to expect that the order of magnitude of the five-loop contribution, from a 
total of 12672 Feynman diagrams [89, 90], is likely to be of 0{a/7rf ~ 7 x IQ-^l 
This magnitude is well beyond the accuracy required to compare with the present 
experimental results on the muon anomaly, but eventually needed for a more precise 
determination of the fine-structure constant a from the electron anomaly [28] . 

7.1.2. Vacuum Polarization Diagrams from Electron Loops 
Vacuum polarization contributions result from the replacement 

■ 9a3 .( QaqA n(^)(g^) 

whereby a free-photon propagator (here in the Feynman gauge) is dressed with the 
renormalized proper photon self-energy induced by a loop with fermion /. Formally, 
with J^Q "(x) the electromagnetic current generated by the /-fermion field at a space- 
time point X, the proper photon self-energy is defined by the correlation function 

t J d'x e*^- (0|T{Ji4) "(x)j£) ^(0)}|0) = - {g-^q^ - q^q^) Il^f\q') .{IQ) 
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In full generality, the renormalized photon propagator involves a summation over all 
fermion loops as indicated in Figure 34 and is given by the expression 

where a is a parameter reflecting the gauge freedom in the free-field propagator (a = 1 
in the Feynman gauge). 



vAAAA/ 



Figure 34. Diagrammatic representation of the full photon propagator in 
Equation (77). 



Since the photon self-energy is transverse in the g-momenta, the replacement in 
Equation (75) is unaffected by the possible gauge dependence of the free-photon 
propagator. 

On the other hand, the on-shell renormalized photon self-energy obeys a dispersion 
relation with a subtraction at = (associated to the on-shell renormalization) ; 
therefore 

q'^ Jo t t — q"^ IT 

where ;^Imn'^-^)(t) denotes the /-spectral function, related to the one-photon e~^e~ 
annihilation cross-section into f'^f~ (see the illustration in Figure 35 ) as follows: 

a(t)e+e-^/+/- = ^^Imn(/)(t). (79) 
More specifically, at the one-loop level in perturbation theory 

-lmU^f\t) = ^i./^TM ( 1 + 2^^] e{t - Am)) . (80) 
Tx 7r3V ^\ ^/ 

The simplest example of this class of Feynman diagrams is the one in Figure 36, 
which is the only contribution of this type at the two-loop level, with the result [91, 81] 



af{m^/m^) 



1) fiVn^-- + 

3/ \2J m, 36 




Pi 

Vacuum polarization contributions from fermions with a mass smaller than that of the 
external line (me <^ m^) are enhanced by QED short-distance logarithms of the ratio 
of the two masses (the muon mass to the electron mass in this case), and are therefore 
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Im 



Figure 35. Diagrammatic relation loctwccn tlic spectral function and the cross-section 
in Equation (79) . 




Figure 36. Vacuum Polarization contribution from a Small Internal Mass 



very important, since In j ~ 5.3. As shown in Reference [92], these contributions 
are governed by a Callan-Symanzik-type equation 

where I3{a) is the QED-function associated with charge renormahzation, and a^(^,Q;) 
denotes the asymptotic contribution to from powers of logarithms of ^ and constant 
terms only. This renormahzation group equation is at the origin of the simplicity of the 
result in Equation (81) for the leading term: the factor 2/3 in front of In ^ comes from 
the first term in the expansion of the QED /3-function in powers of ^ [94]: 




while the factor 1/2 in Equation (81) is the lowest-order coefficient of a/vr in 
Equation (67), which fixes the boundary condition (at = me) to solve the differential 
equation in Equation (82) at the first non-trivial order in perturbation theory, i.e., 
0{^)^. Knowing the QED /3-function at three loops, and a/ (from the universal class 
of diagrams discussed above) also at three loops, allows one to sum analytically the 
leading, next-to-leading, and next-to- next-to-leading powers of In m^/me to all orders 
in perturbation theory [92]. Of course, these logarithms can be re-absorbed in a QED 
running fine-structure coupling at the scale of the muon mass Q;QED('^Ai)- Historically, 
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the first experimental evidence for tlie running of a coupling constant in quantum field 
theory comes precisely from the anomalous magnetic moment of the muon in QED, 
well before QCD and well before the measurement of a{Mz) at LEP[93] (a fact which 
unfortunately has been largely forgotten). 

The exact analytic representation of the vacuum polarization diagram in Figure 36 
in terms of Feynman parameters is rather simple: 

2 



a(^)K/m,)= (^) [' ^{1 - x){2 - x) f dyy{l - y) ^ . (84) 

XTT/ Jo X Jo 1 + !I^-Lz£_-L_ 



This simplicity (rational integrand) is a characteristic feature of the Feynman parametric 
integrals, which makes them rather useful for numerical integration. It has also been 
recently shown [95] that the Feynman parameterization when combined with a Mellin- 
Barnes integral representation is very well suited to obtain asymptotic expansions in 
ratios of mass parameters in a systematic way. 

The integral in Equation (84) was first computed in ref. [96]. A compact form of 
the analytic result is given in ref. [97] , which we reproduce below to illustrate the typical 
structure of an exact analytic expression (p = mg/m^): 



25 1 



36 



y-21npln -Li2(p2) 



5p2 



In p In f 

2 ll + p^ 



Li2(p) + Li2(-p) 



^5) 



Numerically, using the latest CODATA [98] and PDG [99] recommended value for the 
mass ratio m^/rrie = 206.768 2838(54), one finds 

al^\m^/m^) = 1.094 258 3111 (84) (^"^ , (86) 



TT 



where the error in the last two significant figures is the one induced by the present error 
in the mass ratio. 

At the three-loop level, there is the diagram in Figure 37 generated by the insertion 
in the Schwinger diagram in Figure 31 of the vacuum polarization due to two electron 
loops [100], 




Figure 37. Three Loop contribution from Vacuum Polarization corrections due to 
two Electron Loops 
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as well as the three diagrams in Figure 38 generated by the proper fourth-order vacuum 
polarization with one electron loop [101]. 





Figure 38. Three-loop contribution from tlie proper fourtli-ordcr vacuum polarization 
corrections due to an electron loop. 



At the same order, one has to add the 14 Feynman diagrams generated by the 
insertion of the lowest-order vacuum polarization correction due to an electron loop in 
the fourth-order vertex Feynman diagrams in Figure 32. This generates the diagrams 
of Figure 39, which we have collected in five subclasses of gauge- invariant diagrams. 
The history and references of the earlier evaluation of these contributions can be found 
in the review article in ref. [102]. The exact analytic expression from these vacuum 
polarization graphs, for arbitrary values of the masses, was completed in 1993 by Laporta 
and Remiddi [103, 105]. The expression is so lengthy that it is not even reproduced in 
the original papers, where instead, the asymptotic expansion in the mass ratio m^/mg 
is given up to very high order. The first few terms in that expansion are 



+ 



1075 25 



9 



In^ 



216 



— 7r^ + -7r^ln2-3C(3) 





"31 




+ 


27 




11 „ 








-im 


216 


9 



^7r2ln2 + C(3) 



In 



;i/2) 



+- 



nif. 



m,. 



3199 



1080 



-TT 



2 7rMn2 

9 



13 



-TT 



+ 



rap 



m,. 



In^ 



m,, 



(87) 



The leading, next-to-leading, and next-to-next-to leading terms coincide with the earlier 
renormalization group calculation of ref. [92]. Numerically, using the exact analytic 
expression one finds [97] 

3 



,(6), 



ejvp 



1.920 455 130 (33) ( - 
\7r 



while using the asymptotic expression in Equation (87) one gets 

3 



,(6), 



eJvp 



1.920.. 



(89) 



which gives an indication of the size of the terms neglected in Equation (87). 

The contributions at the four-loop level with the insertion of vacuum-polarization 
diagrams from electron loops have been carefully analyzed by Kinoshita and Nio in a 
series of papers, the most recent being Reference [108] where earlier references can also 
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Figure 39. Three-loop contributions from vacuum polarization corrections due to 
an electron loop; the diagrams in each line define a subclass of gauge invariant 
contributions. 

be found. They have classified all the possible Feynman diagrams in three groups of 
gauge-invariant contributions: 

(i) Group I consists of all the possible electron-loop-type vacuum-polarization 
insertions in the photon propagator of the Schwinger graph in Figure 31. All 
together, this group consists of 49 diagrams. Typical examples of Feynman graphs 
of this group are shown in Figure 40. 




Figure 40. Examples of four-loop contributions of group I. 

(ii) Group II consists of the 90 diagrams obtained from all the possible electron-loop- 
type vacuum polarization insertions in the proper two-loop vertex diagrams in 
Figure 32. Examples of Feynman graphs of this group are shown in Figure 41. 
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Figure 41. Examples of four-loop contributfons of group II. 



(iii) Group III consists of the 150 diagrams obtained by the lowest-order electron-loop 
vacuum polarization insertion in all possible photon propagators of the proper sixth- 
order vertex diagrams. Examples of Feynman graphs of this group are shown in 
Figure 42. 




+ 




+ 




+ 




Figure 42. Examples of four-loop contributfons of group III. 



All integrals of these three groups have been evaluated numerically, although some of 
them have also been evaluated using their asymptotic expansion in m^/rrif. [103, 104]. 
The numerical results obtained in Reference [108] are 



ajf)(m^/me)vpi = 16.720 359 (20) 



a 



vr 



af{m^/me)^pii = - 16.674 591 (68 
af(m^/me) vpiii = 10.793 43 (414) (^^J . 



(90) 
(91) 
(92) 



The strong cancellation between the vpl and the vpll contributions, as well as their 
individual size, can be qualitatively understood using the Callan-Symanzik framework 
which describes their respective asymptotic contributions, as discussed in Reference [92] . 
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7.1.3. Vacuum polarization diagrams from tau loops. 

The simplest example of this class of diagrams is the Feynman graph in Figure 43, which 




Figure 43. Vacuum Polarization contribution from a Large Internal Mass 



gives a contribution[101] 



,(4), 



1 

45 



O 



m m, 
log ■ 



(93) 



In full generality, internal heavy masses in the vacuum polarization loops (heavy with 
respect to the external lepton line) decouple; i.e., they give a contribution which vanishes 
in the limit of an infinitely heavy mass. This is the reason why the muon anomaly is more 
sensitive to new Physics than the electron anomaly, roughly by a factor {m^/mef' . This 
is also the reason why the QED theory for the electron anomalous magnetic moment is 
more precise, since the errors induced by the masses of the heavy leptons appear at a 
much smaller level. 

From a structural point of view, the Feynman diagram in Figure 43 has much in 
common with the hadronic vacuum polarization contribution which we shall discuss later 
and therefore it deserves some attention. In particular, one would like to understand how 
the simple leading behavior in Equation (93) arises. For this purpose it is convenient 
to reconsider the equations (75) to (80) for the case where the reference fermion is a 
tauon, i.e., / = r. The contribution to the anomalous magnetic moment of the muon 
from a r-loop vacuum polarization insertion can then be viewed as the convolution of 
the r-spectral function ^lmn^'^^(t) with the contribution to the muon anomaly, induced 
by a fictitious massive photon with a free-propagator: 



-^ 9al3 



(94) 



q J q — t 

The overall contribution to the muon anomaly from Figure 43 can then be written as 



,(4), 



71 Jim% t 71 ^ ' Jo X^ + - 

^ mf. V 



X) 



(95) 



explicitly showing the fact that the integrand is positive and monotonically decreasing 
in the integration region 4m^ < t < oo. (This is why the result in Equation (93) is 
positive.) Notice also that, for large values of t/m?^, the integral over the Feynman 



Muon g-2: Review of Theory and Experiment 



74 



parameter x behaves as |^ + O 



Since the threshold in the r- 



[ml/t) \og{t/ml) 

spectral function is much larger than the muon mass (4m^ ^ "^P? we can approximate 
the function [Figure 43] by its leading asymptotic behavior, which results in 

a^^\mjm^) ^ [^] \ T . ^^^Imn^^^t) = ( ^ ) ^--f , (96) 



.71 J ?>JAml t t IT ^ ' Kir J 3 15 my 

where the remaining integral over the r-spectral function is just the slope of the r- 
vacuum polarization at the origin, which can be easily calculated inserting the explicit 
spectral function in Equation (80). This is how the leading behavior in Equation (93) 
arises. Numerically, using the latest CODATA [98]- and PDG-recommended [99] value 
for the mass ratio m^/rrij. = 5.945 92(97) x 10~^, one finds 



a[f)(m^/m^) 



0.000 078 564 (26) (^^^ 



(97) 



using the simple asymptotic result in Equation (96); while using the exact expression 
in Equation (85) (with p = m^/m^), one gets 

al^\m^/mr) = 0.000 078 064 (25) (-Y , (98) 



showing that the leading asymptotic behavior already reproduces correctly the first two 
non-zero digits of the exact (^^^ coefficient. 

At the three-loop level, the r-vacuum polarization diagrams are those of 
Figs. 37,38,39 with the electron replaced by the tau and hence: rrie ^ m,-. The numerical 
contribution, from the exact analytic calculation of Laporta and Remiddi [103, 105], 
gives the result 

3 



a(f)(m./m^)vp = -0.001 782 33 (48) ( - 



a 



(99) 



To this one has to add the mixed vacuum polarization diagrams in Figure 44, a priori 
enhanced because of the electron loop. 





Figure 44. Mixed vacuum polarization contribution at tlic three-loop level. 



A recent analytic evaluation of the asymptotic contributions up to terms of 

2 22 / 2 ^ \ 

In ^ In '"^ and Oi^^] has been made in Reference [95] using a 

mf, \ mi mi. I L J o 



o 



new technique, which appears to be very powerful, to obtain asymptotic expansions 
of Feynman graphs. The result agrees with the asymptotic terms calculated in 
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Reference [106], using the more complex method of regions developed by Smirnov [107]. 
Numerically, this results in a contribution 

af{m^/m^,m^/mr)^^ = 0.000 527 66 (17) (^^^ . (100) 

At the four-loop level, the vacuum polarization contributions from r-loops that are 
still potentially important are those enhanced by electron vacuum polarization loops. 
They have been calculated numerically by Kinoshita and Nio [108, 109, 110] with the 
result 

a(f)(m^/me,m^/m,),p = -0.046 188 (37) (^^^ . (101) 
7.1.4- Light- by- Light Scattering Diagrams from Lepton Loops. 

It is well known that the light-by-light diagrams in QED, once the full set of gauge- 
invariant combinations is considered, do not require new renormalization counter-terms. 
Because of that, historically, it came as a big surprise to find out that the set of diagrams 
in Figure 45, when the lepton / in the loop is an electron, produced a contribution 




Figure 45. Light-by-light scattering contribution from an Internal Lepton I. 

proportional to log(m^/me); and, in fact, with a large coefficient [111], which was first 
evaluated numerically. Much progress has been made since then, and this contribution 
is now known analytically for arbitrary values of the lepton masses[112]. 

From a theoretical point of view, there are some general basic features of the leading 
light-by-light contributions to the muon anomaly, which will also be useful later for the 
discussion of the hadronic light-by-light contribution, and which we wish to recall here. 
Basically, one is confronted with a vertex function (see Figure 45 for reference to the 
routing of momenta jl): 

tt We use the following conventions: {7^, j^} = Sg^^, with 17^1^ the Minkowski metric tensor of signature 
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2,Pij = ^e 



(27r4) y (27r4) /tfA;^! 



7"(;^2+ ^2 - m^y'Yi^i- h - m,)-^Y , (102) 



where n|j'^^^(g, /ci, /cs, /C2) denotes the off-shell photon-photon scattering amplitude 
(g + A;i + A;2 + ^3 = 0) induced by the lepton /-loop, and q = p2—pi- An ingenious use of 
gauge-invariance, first proposed by the authors of Reference [111], makes the extraction 
of the anomalous magnetic moment much simpler, at the same time explicitly removing 
the spurious logarithmic ultraviolet divergence inherent to each individual photon- 
photon scattering amplitude. Current conservation (or external gauge-invariance) gives 
rise to the Ward identity 

q^Il^ilpA(l,h,h,k2) = 0, (103) 

from which one obtains by differentiation 



Therefore, one can write 



r?(P2,Pi) = g'r2(P2,Pi) 



with 



r 



(2vr4) 



d^ki f d^k2 ^Ilfl^AQ,h,k3,k2) 



l2l2l2 
"-1 "-2 "-3 



(104) 



(105) 



(106) 



The muon anomaly is then extracted from this expression, via the (formally) simple 
projection 

af{mJmiU = ^—iiUT^ + m^)[Y.Y]{^+m,)T%{p,p)] . (107) 

The first few terms in the asymptotic expansion of aff\m^/mi)\^i when the lepton 
in the loop is an electron are 



-2i„ , 2 10_2 



7rMn^ + ^7r"-3C(3) 
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(108) 
(109) 

(110) 



i/2)[jn,jiy], 75 = ^7 7 7 7 , and with the totahy antisymmetric tensor e^y 



fiupa 



chosen so that £0123 = +1- 
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explicitly showing the large coefficient of the In— term ( oTt = 6.58), ffist noticed 
numerically by the authors of Reference [111]. Numerically, using the exact expression 
of Laporta and Remiddi [112], one gets 

(m^/me)ixi = 20.947 924 89 (16) . (Ill) 

The contribution when the lepton in the loop in Fig 45 is heavy decouples, as it 
does for the case of the r. The ffist few terms in the asymptotic expansion are 



19 
16 



:C(3) 



161 ("Jh^^Yx? — 
810 \mr / nil, 



-O 



(112) 
(113) 



Numerically, from the exact expression [112], one obtains a contribution 



af{m^/mr)w\ = 0.002 142 83 (69) 



a 



IT 



(114) 



At the four-loop level there are 180 muon vertex diagrams containing a light- 
by-light-scattering electron-loop sub-diagram with second-order radiative corrections. 
They have been calculated numerically [108] using two different techniques for the more 
complicated subclasses, with the total result 

4 



a^«Hm^/me)w = 121.843 1 (59) - 



a 



(115) 



At the same four-loop level there are still diagrams with mixed electron-loops and tau- 
loops, with one of them of the light-by- light type. Although suppressed by a (m^/m^)^ 
factor, they are still enhanced by a In ^ factor. They have also been calculated 
numerically by Kinoshita and Nio [108], with the result 

4 



a^^\mjme,mjmr)i.y = 0.083 782 (75) 



a 



TT 



(116) 



All together, the purely QED contributions to the muon anomalous magnetic 
moment, including electron- and tau-loops of the vacuum polarization type and/or of the 
light-by-light scattering type, are known to four loops, to an accuracy which is certainly 
good enough for the present comparison between theory and experiment. In Table 8 
we have collected all the results of the QED contributions that we have discussed. The 
numerical results have been obtained using the latest recommended CODATA value for 
the fine-structure constant [98], 

a-^ = 137.035 999 11(46) [3.3 ppb] , (117) 

which, as discussed in the introduction, is dominated by the electron anomalous 
magnetic moment measurements. In fact, the most recent determination of a from 
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Table 8. QED Contributions to the Muon Anomalous Magnetic Moment. The two 
numerical values in the same entry correspond to the CODATA value [98] , and to the 
more recent [28] value of a. The preliminary (^)^ estimate is not listed (see text) but 
it is used in the final comparison of theory with experiment. 



Contribution 



Result in Powers of - 



Numerical Value in 10 Units 



a{2) Eq. (67) 



0.5 ^ 



116 140 973.27 (0.39) 
116 140 972.76 (0.08) 



a|f) Eq. (69) 
af{mjme) Eq. (86) 
a(f)(m^/m,) Eq. (98) 



a|f) (total) 



-0.328 478 965 
1.094 258 311 
0.000 078 064 

0.765 857 410 



00 
08 
26 

27 



{i) 



413 217.62 (0.015) 



ajf) Eq. (70) 



r(f)(m^/me).p Eq. (88) 
mf,/mr)^p Eq. (99) 



,(6) I 



a^^\mf,/me,m^/mr)^p Eq. (100) 



a|fH'^/.Me)ixi Eq. (Ill) 
a(f)(m^/m,)ixi Eq. (Ill) 



ajf) (total) 



1.181 241 46 

1.920 455 13 
-0.001 782 33 

0.000 527 66 
20.947 924 89 

0.002 142 83 

24.050 509 64 



00 
03 
48 
17 
16 
69 

87 



30 141.90 (0.001) 



,(8) 



Eq. (74) 



a(f)(m^/me)vp Eqs. (90), (91), (92) 
a^i^\m^/me,mf,/mr)^pEq. (101) 



, (8) I 



m^/me)\d Eq. (115) 



af^{mf,/me,ni^/mr)\d Eq. (116) 



a|f) (total) 



-1.728 3 


(35 


10.839 2 


(41 


-0.046 2 


(00 


121.843 1 


(59 


0.083 8 


(01 



(?) 



(?) 



130.991 6 (80) (f)' 



381.33 (0.023) 



^(2+4+6+8) (QED) 



116 584 714.12 (0.39)(0.04) 
116 584 713.61 (0.08)(0.04) 



the comparison between QED and a new measurement of the electron (^-factor [20] gives 
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the result [28] 

a^^ = 137.035 999 710(96) [0.70 ppb] , (118) 

which is quite a remarkable improvement. The numerical results using this more precise 
value of a as an input are also given in the Table 8. Use of this new value only affects 
the determination of the lowest-order term a^^-* (the second-line value) and hence the 
total at the bottom in Table 8. 

Notice that the errors in a|,'^^ (total), a|f^ (total) and alf'' (total) are very small and 
dominated by the error in the r-mass, and in the /i-mass to a lesser degree. The first 
error in the total sum a^^"'"^"'"^"'"^^(QED), which is still the dominant error, is the one 
induced by the error in the determination of the fine-structure constant; the second 
error is the one induced by the lepton mass ratios. 

The most recent estimate of the five-loop contribution by Kinoshita [89, 109] gives 
a result 

a^^'\m^/m, , m^/m,),p ^nd ixi ^ 663 (20) (^^^ = 4.48 (0.14) x 10"" . (119) 

Because this estimate is not yet at the level of rigor of the lowest-order determinations, 
we have not included it in Table 8; however, we shall take it into account when making 
a final comparison between theory and experiment. 

The question that naturally arises is whether or not the discrepancy between the 
experimental result on the one hand, and the QED contribution from leptons alone which 
we have discussed, can be understood in terms of the extra hadronic and electroweak 
contributions predicted by the Standard Model. This will be the subject of the following 
subsections. 



7.2. Contributions from Hadronic Vacuum Polarization 

All calculations of the lowest-order hadronic vacuum polarization contribution to the 
muon anomaly (see Figure 46) are based on the spectral representation [113] 

al'\RU = - / --Imn(^^) t / dx \ ^ , 120 

^ 71 Jo t 71 JO X^ + -^(1—X) 

with the hadronic spectral function ^Imn''^^(t) related to the one-photon e^e~ 
annihilation cross-section into hadrons as follows: 

47r2a 1 . 

As already explained in the previous subsection when discussing the contribution from 
vacuum polarization due to a heavy lepton. Equation (120) results from the replacement 

in the free-photon propagator of the lowest-order QED diagram in Figure 31 by the 
one corrected by the proper hadronic photon self-energy in Figure 46. Since the photon 
self-energy is transverse in the g-momenta, the replacement is unaffected by the gauge 



a(t)|e+e-^{7Hhadrons} = ^--^^(t) . (121) 
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Hadrons 

Figure 46. Hadronic Vacuum Polarization Contribution 



dependence in the free-photon propagator. The on-shell renormahzed photon self-energy 
obeys a dispersion relation, with a subtraction at = associated with the on-shell 
charge renormalization; therefore 

Jo tt-q^Ti ^ ' ^ ' 

The fact that only one subtraction is needed in this dispersion relation follows from the 
QCD short- distance behavior of Ii^^\q^) in the deep Euclidean region (— — > oo). The 
r.h.s. of Equation (120) can thus be viewed as the convolution (the integral over t) of 
the hadronic spectral function with the contribution to the muon anomaly, induced by 
a fictitious massive photon {w? = t) with a free-propagator: 

g2 ) q^-t 

This massive photon contribution to the muon anomaly [114], which we shall denote by 
(f) K^^\t/ml), results in a simple Feynman parametric integral (the integral over x in 
Equation (120)): 



The integrand on the r.h.s. of Equation (125) is positive and monotonically decreasing 
in the integration region < t < oo. This is why the lowest-order hadronic vacuum 
polarization contribution to the muon anomaly is positive. This integral was first 
evaluated analytically in Reference [114]. A convenient representation for numerical 
(and analytical) evaluations, which has often been used in the literature, is the one 
given in Reference [116]: 



For t > ^mi, and with = Jl and x 



j^{2) ( _L] = 13,2^2 - x') + i^ + '^^Y^'''^ [ln(l + x)-x+ -x'] + ^x' In x .(126) 
\ mf, / 2 x^ 2 1 — X 



/J., 



Let us also remark that, for large values of t/rnj^, the function K^'^\t/m'j^) decreases as 
ml/t: 



1 rn^ 

lim K<^^\t/ml) = -^ + 



^l/t) log(t/m^ 



(127) 
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and this is why the bulk of the contribution to a^[Fig. 10] comes from the low-energy 
region- in particular from the prominent lowest p-resonance. The shape of the function 
K^'^\t/m?^), as well as its leading asymptotic behavior in Equation (127), are plotted in 
7.2. One can see from this figure that, already at the p-mass value t = M^, the leading 




Figure 47. Behavior of tlic function K^^^t/mj^) in Eq. (125). The dotted blue hne 
represents its leading asymptotic behavior in Eq. (127). 



asymptotic behavior is a rather good approximation to the exact function. In fact, as 
first pointed out in Reference [115], the leading asymptotic behavior of the function 
K^'^\t/m'j^) provides an upper bound to the size of a^[Fig. 10]: 

2 /9^m(g2^^ 
V l- 



,(4) 



(H) 



a 1 



VP < --m- 

TT O 







oo 11 

rft--ImnW(t) 



a 1 



tt3 



-m, 



:i28) 



In other words, to lowest-order in powers of the fine-structure constant, the hadronic 
vacuum-polarization contribution to the muon anomaly has an upper bound which is 
governed by the slope at the origin of the hadronic photon self-energy, a quantity which 
is one of the low-energy constants (a constant of 0{p^)) of the effective chiral Lagrangian 
of QCD. 

From a theoretical point of view, it is more convenient to convert Equation (120) 
into an integral over the hadronic photon self-energy in the Euclidean region, i.e.. 



an integral over 



-q^ > instead of an integral over the Minkowski region 



t > 4m^. This is easily achieved combining Equation (120) and the dispersion relation 
in Equation (78), with the result [116, 102]: 



,(4) 



(H) 



vp 



a 

TV Jo 



dx{l 



X) 



.n(H) 



X 



-m, 



X 



(129) 
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Another useful representation, which follows from this one by partial integration, is the 
one in terms of the Adler-function 

with the result 

= ^\ lo - - "^^'"^ (t^^^O ■ ^'''^ 

These representations are better adapted to theoretical analysis and, eventually, to 
lattice QCD numerical evaluations. In fact, some exploratory lattice QCD work in this 
direction has already started [120]. 

It is possible to make a quick analytic estimate of a|f^(H)vp using a Large- A^'c QCD 
framework (see, e.g. References [117, 118] for review articles and Reference [119]). 
The spectral function in the minimal hadronic approximation (MHA) to Large- A^c QCD 
consists of a lowest narrow state (the p) plus a perturbative QCD continuum starting 
at a threshold Sq, yielding a simple parameterization for the Adler function (for three 
flavors) : 

A^^HQ') - -e' UfM-^ + fl + ^^^^ + ■ ■ (132) 

where jp denotes the p to vacuum coupling, related to the electronic width of the p as 
follows: 

rp-.e+e- = ^/p'M, , (133) 

and the dots stand for higher-order terms in powers of a^. The short-distance scale Sq 
denotes the threshold of the perturbative QCD (pQCD) continuum. The actual value 
of So is fixed from the requirement that, in the chiral limit, there is no 1/Q^-term in 
the short-distance behavior of the Adler function, a property which follows from the 
operator product expansion (OPE). This leads to the constraint 

.9 9 A'c 4 / 3A'cas(so) \ 

2fMl = — ^-so 1 + ^ ^ + ■ ■ ■ • 134 

■'P P 167r2 3 V 8 TT ; ^ ' 

Inserting the large- A"c ansatz given in (132), with = j^^ri^ and the constraint (134) 
in Equation (131), results in the simple estimate 

a(f)(H)^p ~ (57.2 ± 15.0) x 10"^ , (135) 

9 

where we have set = 2-^, as predicted by the MHA to Large-A'c QCD, with Fq the 
pion to vacuum coupling in the chiral limit with: Fo ~ 87 MeV and Mp ~ 750 MeV. 
The error in Equation (135) is a generous estimate of the systematic error of the 
approximations involved. One can, in principle, obtain a much more refined theoretical 
evaluation, but not with the high precision which is nowadays required. Regrettably, 
the present lattice QCD estimates [120] are not more precise than this simple large- A'c 
QCD estimate. 
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The phenomenological evaluation of the hadronic vacuum polarization to the 
muon anomaly, using a{e'^e~ — > hadrons) data, goes back to earlier papers in 
References [121, 37, 122] with the results quoted in the first part of Table 9. (The 
question-mark error in Reference [121] is due to ignorance at that time of the high- 
energy contributions, since this determination was prior to QCD.) The advent of more 
and more precise data on e"*" — e~ annihilations, as well as the new — 2 results from 
the Brookhaven E821 experiment, have motivated more refined re-evaluations of these 
hadronic contributions, which also incorporate QCD asymptotic behavior properties. 
In the second part of Table 9 we collect recent detailed determinations including, in 
particular, the re-analysis of earlier measurements by the CMD-2 detector at the VEPP- 
2M collider in Novosibirsk [123]. It is also possible to measure ^lmll^^\t) in the region 
using hadronic r-decays [124]. Two representative results, which partly use 
r-data, and which also take into account corrections to the conserved-vector-current- 
symmetry (CVC) limit where the spectral functions from r-decays and e^e~ are directly 
related, are quoted in the third part of Table 9. As a test of the CVC hypothesis, the 
e~^e~ data can be used to predict the r branching fractions which, unfortunately, leads 
to branching ratios in disagreement with the data from LEP and CLEO [125]. As 
additional e+e~ and r data have become available, this discrepancy has increased [132]. 

A preliminary update of the hadronic vacuum polarization, contributions reported 
by M. Davier [132] at the Tau06 conference incorporates the new CMD-2 results [129, 
130], new SND results [131], as well as BaBar data for some exclusive channels other 
than vr+vr". Hagiwara, et al., [133] have also included these new data, as well as the 
KLOE data, in a recent reanalysis and obtain a similar result. We quote these recent 
evaluations at the bottom part of Table 9. The overall results quoted in this Table 
illustrate the progress which has been made in this field, especially in the last few years. 
Note that the new data allow for a significant reduction of the experimental error as 
well as the error due to radiative corrections. 

Several remarks, however, concerning the results in Table 9 are in order: 

• There still remains a problematic discrepancy between the most recent e~^e~ 
determinations and those using r-hadronic data. This prevents one from using 
a straightforward average of these different determinations. 

• Various sources of isospin breaking effects, when transforming r-hadronic spectral 
functions to the required e~^e~ spectral function, have been identified. A 
straightforward correction is the one due to the electroweak short-distance 
contribution to the effective four-fermion coupling r Vrdu [134, 135, 136]. 
Effects due to the mass differences between the ti^ and the 7r°, as well as mass 
and width differences between the charged and the neutral p, have also been 
considered [124, 137]. 

• A more problematic issue is the one of radiative corrections for a specific channel. 
Calculations concerning the mode z/,-7r~7r° have been reported in References [138, 
139]. 
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Table 9. Compilation of Results on a/i (H) 



Authors 


Contribution to x 10^^ 


Gourdin, de Rafael [121] (1969) 
Bailey et al. [37] (1975) 
Calmet-Narison-Perrottet-de Rafael [122] (1976) 


6500 ± 500 ± (?) 
7300 ± 1000 
7020 ± 800 


Davier et al, [e+e^-data] [125] (2003) 

Jegerlehner [126] (2004) 
de Troconiz, Yndurain [127] (2004) 
Hagiwara et al, [128] (2004) 


6963 ± 62,.p ± 36,.d 

6948 ± 86 
6944 ± 48,.p ± 10,.d 
6924 ± 59..P ± 24.,, 


Davier et al, [r-data] [125] (2003) 
de Troconiz, Yndurain [127] (2004) 


7110±50,.p±8,.d±28sc.(2) 
7027 ±47.., ±10.., 


Davier, [e+e'-data] [132] (2006) 
Hagiwara et al, [e"'"e~-data] [133] (2006) 


6908 ± 30,,p ± 19.,, ± 7qcd 
6894 ± 42,,, ± 18.,, 



• The most precise e"'"e~ -k^-k^ data come from the CMD-2 detector, which are 
now in agreement with data from SND, but not with the recent data from KLOE 
at Frascati. In fact, the two most recent compilations[132, 133] differ on the use of 
KLOE's results, which agree with CMD on the integral but not the shape of the 
pion form factor. The CMD-2 results have been corrected for leptonic and hadronic 
vacuum polarization effects, so that the measured final state corresponds to tt+tt", 
including pion-radiated photons and virtual final-state QED effects. 

• By contrast, hadronic r-data are available from the four LEP detectors, CLEO, 
OPAL, L3, and ALEPH, and show good consistency when measuring the r~ — > 
Vt-ti'ti^ branching ratio. 

• One can hope that the forthcoming results from the KLOE detector [140] and 
BaBar [141], using the radiative return technique, will eventually provide the 
necessary consistency check of the e+e~ input. The BELLE r-data, with one 
hundred times the statistics of the LEP experiments, should also provide a test 
of the LEP and CLEO results. 

• Finally, let us comment on the fact that the contribution from the high-energy 
tail to a|f)(H)vp in Equation (120) is well under control because of the asymptotic 
freedom property of QCD. The perturbative QCD prediction is known to next- 
to-next-to-leading order with second-order quark mass corrections included. (See 
Reference [142] and references therein for earlier calculations.) 

7.2.1. Higher- Order Hadronic Vacuum Polarization Contributions 
There are three classes of hadronic vacuum polarization contributions of an overall 
O 



which were first considered in Reference [122]. A first class is the set of 
Feynman diagrams in Figure 48. They are obtained by inserting one hadronic vacuum 
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polarization correction into each of the virtual photons of the purely QED fourth-order 
diagrams in Figure 32. The diagrams in each line in Figure 48 are well-defined gauge- 









Figure 48. Hadronic Vacuum Polarization Corrections to the two- loop graphs in 
Figure 32. 



invariant subsets. One can easily see that their overall contribution, by analogy to our 
discussion of the lowest-order hadronic vacuum polarization, will be governed by an 
integral 



,(6) 



(H)vp 



dt 1 



t 



Imn(H)(t)K(^) ( - , 

t TT V^^. 



(136) 



where K^'^^ contribution to the muon anomaly from a fictitious 

massive photon in the presence of O (^^^ QED corrections. The exact form of the 

function K^^'^ which is a rather complicated expression, was obtained by Barbieri 

and Remiddi [143] and has not been fully checked since then. We reproduce in Figure 49 
the shape of this function and compare it with its leading asymptotic behavior, which 
agrees with the one that has been extracted [144] from an independent effective field- 
theory calculation of anomalous dimensions of composite operators in Reference [145], 



hmir(^)f^U-^lilnaUof< 



(137) 
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0.4 




Figure 49. Behavior of tlie function K^^^t/mj^) in Eq. (136). The dotted blue hue 
represents its leading asymptotic behavior in Eq. (137). 



The fact that the function K^^^ ( ) is negative for 4yU^ < t < oo imphes that, contrary 



to the lowest-order case, the contribution to a|f^(H)vp must be negative. 

A second class of potentially important contributions are the mixed vacuum 
polarization insertions in Figure 50 due to an electron loop and a hadronic loop. [The 
similar contribution where the electron loop is replaced by a r loop decouples as m^/m^ 
and is negligible at the present required level of accuracy.] This contribution, can be 




Figure 50. Second Class of Higher-Order Hadronic Vacuum Polarization 
Contributions. 



expressed double convolution [122]: 



.{-2)f(-i)log 



where, in underbraces, we have indicated the leading asymptotic behaviors in rnj^/t and 
ml/m'jj^- indicating that this contribution must be positive. The question remains: which 
of the two contributions from Figure 48 (which, we concluded, is negative) and from 
Figure 50 (which is positive) dominates? A simple answer to this can be found, without 
doing the full calculation [146], by combining the asymptotic behaviors in Eqs. (137) 
and (138): 
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vp 



- / ^-Imn t 

71 J Jams t t tt 



23, t 2 ml 
In — + - In 



ml 2 IQtx'^F, 



TT 



2t?2 




Af 2 3 Ml 



23, M2 

Ill TT 



m. 



9 mf 



(139) 
(140) 



where in the second hne we have use the QCD Large- A^c dominance of the p contribution. 
The overall result is indeed negative, because of the larger weight in the hadronic term. 
This cancellation, however, also shows the need of a more refined estimate. Indeed, the 
most recent evaluation, using the same e"'"e~ data as in the lowest-order calculation, 
gives [147] 

m,. 



,(6) 



H and — , H 



mp 



-97.9 ± 0.9_ ± 0.3„.H X 10 



-11 



vp 



(141) 

101 ± 1 X 10-" if the r 



there are still two extra 



This determination goes down to a|f^(H and m^/me , H)vp 
data mentioned previously are used. 

As first mentioned in Reference [122], at O 
contributions from hadronic vacuum polarization which have to be considered. They are 
shown in Figure 51. The evaluation of the diagram in Figure 51(A) can, in principle, 
be achieved by lumping together in the lowest-order hadronic vacuum polarization, 
the effect of replacing the one-photon exchange cross-section in Equation (121) by 
the physical cross-section which incorporates the full hadronic vacuum polarization 
(i.e., not just the proper hadronic self-energy). The contribution from the diagram 
in Figure 51(B) can also be incorporated in the lowest-order evaluation, if the one- 
photon-exchange hadronic cross-section also takes into account final states with one- 
photon emission as well as the one photon radiative corrections to the hadronic vertex. 
Unfortunately, extracting the latter from the full data, involves a certain amount 
of model dependence, with uncertainties which are difficult to ascertain from the 
discussions in the literature. In principle, these uncertainties should be reflected in 
the radiative correction error quoted in Eq. (141). 





(B) 



Figure 51. Further contributions from higher-order hadronic vacuum polarization. 
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7.3. Contributions from Hadronic Light-by-Light Scattering 

These are the contributions illustrated by the diagrams in Figure 52. All the estimates of 
these contributions made so far are model dependent. There has been progress, however, 
in identifying the dominant regions of virtual momenta, and in using models which 
incorporate some of the required features of the underlying QCD theory. The combined 
frameworks of QCD in the l/A'^c-expansion and of chiral perturbation theory[148] have 
been very useful in providing a guideline to these estimates. 



X 



1^ 




+ Permutations 



Figure 52. Hadronic Light-by-Light Contributions 

So far, the only rigorous result in this domain has come from the observation [149] 
that, in large-A^c QCD and to leading-order in the chiral expansion, the dominant 
contribution to the muon g — 2 from hadronic light-by- light scattering comes from the 
contribution of the diagrams which are one-particle (Goldstone-like) reducible. The 
diagrams in Figure 53(a) produce a In^ (/x/m)-term modulated by a coefficient which 
has been calculated analytically [149], with the result: 

Here, denotes the pion coupling constant in the chiral limit (F^ ~ 90 MeV); the 
/i-scale in the logarithm is an arbitrary ultraviolet (UV)-scale, and the m-scale is an 
infrared (IR) mass, either or if only the leading term in Equation (142) is 
known. 




(a) (b) (c) 

Figure 53. One Goldstone Reducible Diagrams in Chiral Perturbation Theory 
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The arbitrariness on the UV and IR scales in Equation (142) would be removed, if 
one knew the terms linear in log/i from Figure 53(b), as well as the constant K(/i) 
from the local counter-terms generated by Figure 53(c). Unfortunately, neither the 
determination of the coefficient of the log/x-term, nor the constant K(/i)-term, can be 
made in a completely model- independent wayff. 

The nice feature about Equation (142) is that, as discussed in Reference [149], it 
plays a fundamental role in fixing the overall sign of the hadronic light-by-light scattering 
contribution to the muon g—2. Indeed, in the various hadronic model calculations of this 
contribution, there appear hadronic scales (like the p-mass), which act as a UV- regulator 
and play the role of /i in Equation (142). Therefore, letting the hadronic leading mass- 
scale become large, and provided that the model incorporates correctly the basic chiral 
properties of the underlying QCD theory, must reproduce the characteristic universal 
ln^(/i) behavior with the same coefficient as in Equation (142). This test, when applied 
to the most recent earlier calculations [152, 151] (prior to the Knecht-Nyffeler calculation 
in Reference [153]), failed to reproduce the sign of the coefficient of the ln^(/i)-term in 
Equation (142). The authors of References [152, 151] have later found errors in their 
calculations which, when corrected, reproduce the effective field theory test. Their 
results now agree with the Knecht-Nyffeler calculation [153] on which we next report. 

The Knecht-Nyffeler Calculation [153] 

In full generality, the pion pole contribution to the muon anomaly has an 
hadronic structure, as represented by the shaded blobs in Figure 54. The authors 
of Reference [153] have shown that, for a large class of off-shell 71^77 form 
factors (which includes the \aige-Nc QCD class), this contribution has an integral 
representation over two Euclidean invariants and Q2 associated with the two 
loops in Figure 54: 

roo poo 

a^[Fig. 54]= / dQl dQl W{QlQl)n{QlQl) , (143) 
Jo Jo 

where W(Q^, Ql) is a skeleton kernel which they calculate explicitly, and 'H{Q\, Ql) 
is a convolution of two generic ^TrO-y*^* (^i, 5'|) form factors. In Large- iVc QCD, 

with the sum extended to an infinite set of narrow states. In practice, the calculation 
in[153] has been done by restricting the sum in Equation (144) to one- and two- 
vector states, and fixing the polynomial Cij(g^,g|) from short-distance and long- 
distance QCD properties. This way, they obtained the result 

ajf)(7r°)M = (5.8 ± 1.0) x 10"^° , (145) 

where the error also includes an estimate of the hadronic approximation. Further 
inclusion of the rj and 77/ states results in a final estimate 

a(f)(7r° + 7] + T^Oixi = (8.3 ± 1.2) x 10"^° , (146) 

ft Reference [149] provides a discussion of this point using a rcnormalization group approach. EsscntiaUy 
the same arguments have been subsequently emphasized in Reference [150]. 
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+ Permutations 



Figure 54. Hadronic light-by-light from a 7r° intermediate state. 

in agreement with the earher calculations by Bijnens at al [151] and Kinoshita et 
al [152] (after correcting for the overall sign). 

Comment on the Constituent Quark Model ( CQM) 




+ Permutations 



Figure 55. Hadronic light-by-light in the constituent quark mode 



We would like to comment on an argument which is often used in favor of the 
Constituent Quark Model as a "simple" way to fix the sign and size of the hadronic 
light-by-light scattering contribution to the muon g — 2. Since the argument has 
even been advocated in some publications, we feel obliged to refute it here, with the 
hope that it will stop further confusion. The constituent quark model contribution 
from the diagram in Figure 55, can be easily extracted from the work of Laporta 
and Remiddi [112] for a heavy lepton contribution (see Equation (112) above), with 
the result 
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0.616 



Seen from a low-energy-effective-field-theory point of view, the constituent quark 
mass Mq in the CQM should provide the UV-regulating scale. However, the model 
is not a good effective theory of QCD and, therefore, it fails to reproduce the 
characteristic QCD In^ Mq behavior when Mq is allowed to become arbitrarily 
large; the CQM contribution decouples in the large Mg-limit. Incidentally, the 
contribution from free u and d quarks with a small mass m <^ /i goes as 
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~ Nc^f X l^^l'^^' which is also incompatible with the QCD result in 
Equation (142). Therefore, arguments based on the fact that the CQM (and/or 
pQCD) gives a positive contribution are certainly "simple, " but also incorrect. 
Notice however, that, contrary to the naive CQM, the constituent chiral quark 
model of Georgi and Manohar[154] (see also Reference [155]) does indeed reproduce 
the correct In^ Mq behavior in the Mq — > oo limit. This is because, in these models, 
the Goldstone particles couple with the constituent quarks in a way which respects 
chiral symmetry, and the pion pole diagram appears then explicitly. The same 
happens in the extended version of the Nambu-Jona-Lasinio model [156]. These 
models, however, suffer from other diseases [157], and therefore they are not fully 
reliable to compute the hadronic light-by-light scattering contribution. 

Hopefully, hadronic models of the light-by-light scattering contribution which 
respect the QCD constraint in Equation (142) will be progressively improved, so as to 
incorporate further and further QCD features; in particular, short-distance constraints, 
following the lines discussed in References [117, 118, 158]. An interesting contribution in 
this direction has been reported in Reference [159]. Unfortunately, as recently discussed 
in Reference [160] their numerical evaluation is incomplete and model dependent with 
largely underestimated errors. We believe that, at present, one can only claim to know 
the hadronic light-by-light scattering contribution with a cautious error, which takes into 
account the uncertainties from the integration regions which remain model dependent. 
While awaiting further improvement, the educated value one can quote at present, based 
on the combined work of References [153] and[152, 151] (appropriately corrected) as well 
as Reference [159], is 

a(f)(H)M = (ll±4)xlO-i°. (148) 
7.4- Electroweak Contributions 

The leading contribution to the anomalous magnetic moment of the muon from the 
electroweak Lagrangian of the Standard Model, originates at the one-loop level. The 
relevant Feynman diagrams, which for simplicity we draw in the unitary gauge, are 
shown in Figure 56, where we also indicate the size of their respective contributions. 
The analytic evaluation of the overall contribution gives the result [22] 




+ 38.9x10 -19.4x10 <3.3x10 " 



Figure 56. Weak interactions at the one-loop level 
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19.48 X 10 
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(149) 



where the weak mixing ang le is defined by sin^ Ow = I - M^/M^ = 0.223 (for a Higgs 
mass Mh = 150 GeV), and where Gp = 1.16639(1) x 10~^ is the Fermi constant. Notice 
that the contribution from the Higgs boson, shown in parametric form in the second 
line, is of O In U which is already very small for the present lower bound 

of Mh- Closed analytic expressions of the electroweak one-loop contribution can be 
found in Reference [161] and references therein. 

Let us recall that the present world average experimental error in the determination 
of the muon anomaly is[26] ^a^\^^^= ±6.3 x 10~^° , and, hoping for a continuation of 
the BNL experiment, it is expected to be further reduced. The a priori possibihty that 
the two-loop electroweak corrections may bring in enhancement factors due to large 
logarithmic factors, like ln(M|/m^) ^ 13.5, has motivated a thorough theoretical effort 
for their evaluation, which has been quite a remarkable achievement, and which we next 
discuss. 

It is convenient to separate the two-loop electroweak contributions into two sets 
of Feynman graphs: those containing closed fermion loops, which we denote by 
a^'^^2''(fcrm), and the others which we denote by a^^^'^\ho&). With this notation, the 
electroweak contribution to the muon anomalous magnetic moment becomes 



-.W 



a"" ^-^(bos) 



a"" ^"^(ferm) 



(150) 



with a); the one-loop contribution in Equation (149). 



Bosonic Contributions to aj^^^'^\hos). 

The leading logarithmic terms of the two-loop electroweak bosonic corrections have 
been extracted using asymptotic expansion techniques. In fact, these contributions 
have now been evaluated analytically, in a systematic expansion in powers of 



sin'^ 6w, up to 0[{sm^ Ow)"^] , where log terms, log terms, ^^log 



Mi 



Ml 



terms, terms, and constant terms are kept [162]. Using siv? 6w = 0.224 and 
Mh = 250 GeV, the authors of Reference [162] find 
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(151) 
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Fermionic Contributions to a^^*^^^ (ferm) . 

The discussion of the two- loop electroweak fermionic corrections is more dehcate. 
Because of the U{1) anomaly cancellation between lepton loops and quark loops 
in the electroweak theory, one cannot separate hadronic from leptonic effects 
any longer in diagrams like the ones shown in Figure 57, where a triangle with 
two vector currents and an axial-vector current appears (the so-called VVA- 
triangle). Individually, the lepton- loop and quark-loop contributions are each 
gauge-dependent; depending on the gauge choice, they can even lead to UV- 
divergent contributions. Only the sum of contributions within each family of leptons 
and quarks is free from these ambiguities. As first discussed in References [163, 164], 
it is this anomaly cancellation which eliminates some of the large logarithms that 
were incorrectly kept in a previous calculation in Reference [165]. It is therefore 
appropriate to separate the two-loop electroweak fermionic corrections into two 
classes. One is the class arising from Feynman diagrams like the ones in Figure 57, 
where a subgraph with a VVA-triangle of leptons and quarks appears, including the 
graphs where the Z lines are replaced by lines, if the calculation is done in the 
^z-gauge. We denote this class by aj^^^'^\l,q) . The other class is defined by the 
rest of the diagrams, where quark loops and lepton loops can be treated separately, 
which we call a^^*-^^ (ferm-rest) i.e., 

aj^(2)(fcrm) = aj^(2)(/, q) + af^ ^^\icru.-,csl) . 







1 






< ^ 


z! 








1 


1^ 




Figure 57. Two-loop electroweak diagrams generated by the 77Z-Triangle. There 
are siniilar diagrams corresponding to the /i ; c, s and t ]t,b generations. 



The contribution from a^^''^'*(ferm-rest) brings in m'l/M'^ factors. It has been 
estimated, to a very good approximation, in Reference [164] with the result, 

Gf ml a 



^(^)(ferm-rest) = ^ ^ - X (-21 ± 4) , (152) 
V2 Svr^ vr 



where the error here is the one induced by diagrams with Higgs propagators with 
an allowed Higgs mass in the range 114 GeV < Mh < 250 GeV. 
Concerning the contributions to a^"^ ^ q) ^ it is convenient to treat the three 
generations separately. The contribution from the third generation can be 
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calculated in a straightforward way, because all the fermion masses in the triangle 
loop are large with respect to the muon mass, with the result [163, 164] 



aj^(^)(r,t,6) 




_31o M_io M_^io ^ 

Gp ml a 



^ - X (-30.6) . (153) 



The leading terms can be easily obtained using effective field theory tech- 
niques [163]; but terms of O ^^log-^^ and O ^-^^ have also been calculated 

in Reference [164]. There are in principle QCD perturbative corrections to this 
estimate which have not been calculated, but the result in Equation (153) is good 
enough for the accuracy required at present. 

As emphasized in Reference [163], an appropriate QCD calculation when the 
quark in the loop of Figure 1 is a light quark should take into account the 
dominant effects of spontaneous chiral-symmetry breaking. Since this involves 
the u and d quarks, as well as the second-generation s quark, it is convenient 
to lump together the contributions from the first and second generations. An 
evaluation of these contributions, which incorporates the QCD long-distance chiral 
realization [163, 166] as well as perturbative [145] and non-perturbative [166] short- 
distance constraints, gives the result 

aj^(2)(e, /X, rf, c) = % ^ - X (-24.6 ± 1.8) . (154) 



From the theoretical point of view, this calculation has revealed surprising 
properties concerning the non- anomalous component of the VVA-triangle [167], 
with a new set of non-renormalization theorems in perturbation theory [167, 168] 
(see also References [169, 170, 171]). The physical meaning of these perturbation 
theory results when considered in the full QCD vacuum is, however, still 
unclear [168]. 

Putting together the numerical results in Eqs.(151), (152), (153) with the more 
recent result in Equation (154), one finally obtains the value 
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(156) 



where the first error is essentially due to the Higgs mass uncertainty, while the second 
comes from higher-order hadronic effects in the VVA loop evaluation. The overall result 
shows indeed that the two-loop correction represents a sizable reduction of the one-loop 
result by an amount of 21% . This result has prompted the evaluation of the electroweak 
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leading-log three-loop effects of O 



using renormalization group 



Gf ^ /^Vln Mz. 

^ 87r2 U/ rn^ 

arguments similar to the ones discussed in §(7.1.2). The results [172, 145] show that 
those effects are negligible [~ (9(10~^^)] for the accuracy needed at present. 

7.5. Summary of the Standard Model Contributions to the Muon Anomaly 

The present status in the evaluation of the anomalous magnetic moment of the muon 
in the Standard Model can be summarized as follows: 

• Leptonic QED Contributions 

a^(QED) = (116 584 718.09 ± O.Msioops ± 0.08„ ± 0.04„s) x 10^" . (157) 

This is the number one obtains by adding to the final result at the bottom of 
Table 8 (using recent [28] value of a) to the estimate of the five-loop contribution 
in Equation (119). 

• Hadronic Contributions 

— Hadronic Vacuum Polarization with EM-Data 

From the two most recent determinations in Table 9, which take into account 
the new data on e"'"e~ annihilation into hadrons, we get 

a^[HVP(06)] = (6901 ± 42,,, ± 19.., ± 7qcd) x lO^^^ , (158) 
where we have averaged the central values, but kept the largest uncertainties in 
the two determinations. Notice that this result is consistent with the average of 
the central values of the earlier determinations from e"'"e~ data in the second 
part of Table 9, which when keeping the largest error from experiment and 
from radiative corrections give: a^[HVP] = (6945 ± 62,,, ± 36rad) x 10^^^. 

— Higher-Order Hadronic Vacuum Polarization with EM-Data 

a^[HVP h.o.] = (-97.9 ± 0.9,,, ± 0.3,,^) x 10"^^ . (159) 

— Hadronic Light- by- Light Scattering 

a^[HLLS] = (110 ± 40) x 10"^^ . (160) 

• Electroweak Contributions 

a^[EW] = {154:± 2 ±1) X IQ-^^ (161) 

The sum of these contributions, using the HVP06 result in Equation (158) and adding 
experimental and theoretical errors in quadrature, gives then a total 

These determinations are to be compared to the experimental world average in 
Equation (21) 

a"^" = 11 659 2080 (63) x 10"^^ . (163) 

Therefore, we conclude that, with the input for the Standard-Model contributions 
discussed above, one finds at present a 3.4 a discrepancy. 
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8. Future Experimental Possibilities 

The discrepancy between the experiment and the theory gives strong motivation to 
improve the experimental number. An upgraded experiment, E969[173], with goals 
of 

'^syst — 0.14 ppm and (Jstat — 0.20 ppm received scientific approval from from 
Brookhaven Laboratory. Despite the 3.4 o" difference between theory and experiment, 
the new experiment remains unfunded. Such an experiment would stimulate additional 
theoretical work which would further decrease the theoretical error. The ultimate 
experimental sensitivity of a future experiment depends both on the available beam 
intensity, and on the ability to further reduce the systematic errors. 

The systematic and statistical errors by year for E821, along with the goals for E969, 
are given in Table 10. There is a clear improvement in the systematic errors with time, 
and the improvement required for E969 is modest and should be readily achievable with 
planned upgrades in the detector system, electronics, and the field mapping system. In 
order to achieve the required number of muons for E969, a plan to improve the beamline 
in order to increase the stored muon fiux by a factor of five has been proposed. 

Table 10. Systematic and statistical errors in ppm for each of the e821 running 
periods. 



Systematic uncertainty 


1998 


1999 


2000 


2001 


E969 Goal 


Magnetic Field uip (ppm) 


0.5 


0.4 


0.24 


0.17 


0.1 


Anomalous Precession uja (ppni) 


0.8 


0.3 


0.31 


0.21 


0.1 


Statistical Uncertainty (ppm) 


4.9 


1.3 


0.62 


0.66 


0.20 


Total Uncertainty (ppm) 


5.0 


1.3 


0.73 


0.72 


0.25 



The E821 experimental approach, which combines the "magic 7," electrostatic 
focusing, and a fast muon kicker has the potential to push the experimental error below 
the 0.1 ppm level, perhaps as far as about 0.05 ppm before too many technical barriers 
are reached. To improve the statistical error on uJa the muon fiux would have to be 
increased substantially, and presently such a source does not exist. As the muon fiux 
is increased, systematic errors from pulse pile-up and detector instabilities generated 
by high rates become more severe. To improve the Up measurement, additional work 
would have to be done on the absolute calibration, which at present depends on the 
measurement of Phillips, et al.,[65] to relate the magnetic moment of a proton in a 
spherical water sample to that of the free proton. An alternate calibration based 
on a ^He probe has been suggested[174]. To go beyond 0.1 ppm, it would also be 
necessary to further refine the measurement of A+ with a new muonium hyperfine 
splitting measurement [174]. 

Progressing beyond 0.05 ppm might require a new technique. For example, 
increasing 7r by going to a higher muon momentum and increasing fa by increasing 
the magnetic field would reduce the statistical error without requiring more muons (see 
Equation 11). However, a higher momentum would require abandoning the advantages 



Muon g-2: Review of Theory and Experiment 



97 



of the magic gamma. Increasing the magnetic field would require the elimination of the 
iron yoke and pole tips in favor of a magnetic field driven entirely by superconducting 
coils. A proposal to use sector focusing dipoles has been made [175], but the need to 
know J B ■ di for the muons to a precision of 10 ppb with magnetic gradients present is 
challenging, to say the least. 

9. Summary and Conclusions 

Measurement of the muon's anomalous magnetic moment has a history stretching back 
to the late 1950s. In the paper from Nevis Laboratory that presented evidence for parity 
violation in muon decay [30], it was reported that the data were consistent with g = 2. 
In a subsequent paper from Nevis [33], it was demonstrated that the muon anomaly was 
consistent with the QED calculation of Schwinger of ~ a/27i, and thus provided the 
first evidence, confirmed by all existing subsequent data, that the muon behaves like a 
heavy electron. 

The result from E821 at the Brookhaven National Laboratory is sensitive to 
extremely refined effects predicted by the Standard Model, such as the hadronic light- 
by-light scattering and the electroweak contributions. If these two effects were removed 
from the theoretical prediction obtained in Equation 162, the discrepancy with the 
experimental result would be 7.1 cr. If only the electroweak contribution is removed, 
the discrepancy is 5.1 a, clearly showing the sensitivity of the E821 experiment to the 
electroweak contribution. The largest uncertainties in the theoretical prediction of 
come, from the systematic errors in the experimental data from e"*" + e~ annihilation to 
hadrons (in spite of the fact that the accuracy in the present determination (< 1%) 
is quite remarkable); and from the theoretical error in the hadronic light-by-light 
contribution. 

There is progress to be expected as additional e~^e~ data become available. The 
good news is that the vrvr data from two independent experiments at the Novosibirsk 
collider, CMD-2 and SND, are in agreement. There seems to be some disagreement 
between the shape of the pion form factor obtained by the KLOE experiment at Frascati 
and that from the Novosibirsk experiments. Further experimental input in the near 
future is expected on e~^e~ tt+tt" from a new KLOE analysis with a tagged photon 
in the detector, and from BaBar using initial-state radiation. At present the data from 
Novosibirsk dominate the determination of a^[HVP] in the p-meson region, making an 
independent confirmation of the tt tt cross section in this energy region clearly desirable. 

On the other hand, the puzzle centered on hadronic r-decay data versus e'^e~ 
data remains in spite of efforts to quantify the CVC corrections. Problems exist in 
comparisons between the two that are quite independent of the contribution to a^. For 
example the failure to predict the measured r branching fractions from the e^e~ data 
and the conserved vector current hypotheses, at present invalidate attempts to use these 
r data to obtain the hadronic contribution [132]. Work on these issues in r physics is 
continuing. We conclude that until the issues with r-decay data are resolved, the e~^e~ 
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data are preferred to evaluate the lowest-order hadronic contribution to a^. 

One should stress the impressive achievement reached in the calculation of the QED 
contribution, where the dominant error at present comes from the five-loop contribution. 
Thanks to the dedicated efforts of T. Kinoshita and collaborators, this uncertainty is 
likely to go down significantly in the near future. Great progress has also been made in 
understanding the electroweak contribution, which is now fully calculated at the two- 
loop level with an accuracy of 1.5%, an uncertainty dominated by the unknown Higgs 
mass. The main effort should now be concentrated on a more accurate determination of 
the hadronic light-by-light contribution. A goal of 15% accuracy in its determination 
seems possible with a dedicated effort from theorists. 

An improvement in the experimental value of the muon anomaly by at least a factor 
of two appears to be quite feasible, and is the goal of the proposed E969 experiment at 
Brookhaven National Laboratory. Given the current 3.4-standard-deviation difference 
between theory and experiment, and the prospects for the uncertainty in the theoretical 
value to decrease over the next few years, such a measurement appears timely and highly 
desirable. 

Historically the muon anomaly has provided an important constraint on speculative 
theories of physics beyond the Standard Model [23, 24, 25]. As particle physics moves 
into the Large Hadron Collider era in the next few years at CERN, the muon anomaly 
will provide an independent constraint on the interpretation of the discoveries made at 
the LHC. 
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